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SECTION  I 


INTRODUCTION 


SCOPE 


The  JMEM/SS  Model  for  the  Direct  Fire  Computer  Program  was  developed 
in  1973  by  the  Material  Systems  Analysis  Division,  Rock  Island  Arsenal, 
Ill.,  in  conjunction  with  a  JMEM  Methodology  and  Evaluation  Working  Group. 
Its  purpose  is  to  evaluate  direct  fired  projectile  warhead  kill  mechanisms 
against  typical  battlefield  tactical  elements.  The  projectile  warheads 
most  frequently  used  in  the  simulation  model  are  characterized  by  kinetic 
energy  and  shaped  charge  rounds. 

Program  output  includes  mean  target  Pk  data,  mean  time  Pk  data, 
mean  target  Ph  data,  mean  time  Ph  data,  and  the  average  number  of  rounds 
fired  per  case. 

These  results  are  accrued  from  9,604  Monte  Carlo  trials  which  yield 
target  Pk  and  time  Pk  results  for  four  selected  target  kill  criteria 
such  as  mobility  kills,  firepower  kills,  mobility/firepower  kills,  and 
catastrophic  kills.  Personnel  casualties  for  an  APC  (Armored  Personnel 
Carrier)  type  target  can  also  be  evaluated.  The  computational  sequence 
for  a  Monte  Carlo  trial  is  shown  in  Figure  1-1. 


Figure  1-1.  Basic  Trial  Methodology 

The  four  target  kill  criteria  are  associated  with  each  point  target 
for  a  given  input  data  set.  The  target  area  is  described  by  grid  cells 
on  a  target-enveloping  rectangle  with  dimensions  determined  by  the 
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extremities  of  the  target  form  as  shown  in  Figure  1-2.  Input  target  Pk 
coordinates  are  converted  into  the  identifying  kill  (IK)  array  subscripts 
which  correspond  with  the  subscripts  of  the  target-enveloping  grid  area 
arrays  (AX  and  AY).  The  input  provides  a  Pk  map  which  correlates  the 
target  vehicle  grid  cells  with  associated  vulnerability  data.  Thus, 
only  those  grid  cells  which  have  correlating  vulnerability  data  identify 
the  target  form. 


Figure  1-2.  Target  Description 


The  Monte  Carlo  sampling  technique  establishes  locations  for  round 
aimpoints  and  impact  points  for  both  main  rounds  and,  if  applicable, 
spotting  rifle  rounds  (used  as  a  delivery  assist  device).  If  the  target 
is  hit  by  a  main  round,  tests  are  performed  to  ascertain  the  level  of 
damage  inflicted.  Monte  Carlo  sampling  is  also  used  to  determine  firing 
time  factors  contributing  to  the  evaluation  of  time-to-kill . 

The  successive  Pk  and  Ph  trial  results  are  accumulated  and  reduced 
to  mean  values  for  program  output  after  completion  of  all  the  Monte 
Carlo  trials. 
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COMPUTER  REQUIREMENTS 


The  program  accepts  card  input  data,  contains  about  400  source 
statements,  and  requires  approximately  200,000  bytes  of  core.  An 
average  target/weapon  simulation  consumes  about  6  minutes  of  IBM  360 
running  time.  JMEM  groups  constitute  the  major  users  of  this  program 
in  the  field. 

Preparation  of  complete  vulnerability  data  for  a  target  may  prove 
to  be  a  laborious  task  since  the  data  is  not  of  a  general  nature  but 
describes  only  a  specific  attack  elevation  and  azimuth  angle  for  a 
specific  kill  criterion  and  weapon.  Users  can  obtain  the  required 
vulnerability  data  from  BRL,  Aberdeen  Proving  Grounds,  Md. 
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SECTION  II 


MATHEMATICAL  MODEL 


DELIVERY  ERROR  MODEL 

A  Monte  Carlo  sampling  technique  is  used  to  estimate  the  projectile 
aim  and  impact  point  locations.  It  is  assumed  that  the  true  aimpoint  of 
the  projectile  is  normally  distributed  about  the  intended  aimpoint,  and 
the  projectile  impact  point  is  normally  distributed  about  the  true 
aimpoint.  ' 

Main  Round 

The  aimpoint  (center  of  impact  or  Cl)  of  the  first  main  round  for 
each  trial  is  distributed  probabilistically  about  the  target.  The  aim- 
point  of  succeeding  rounds  fired  within  a  sample  either  remains  the  same 
or  is  modified  by  a  sensed  miss.  Thus,  during  the  course  of  the  trials, 
the  group  of  impact  points  will  tend  to  be  distributed  about  the  target 
according  to  a  probability  distribution  which  is  dependent  on  the  aiming 
accuracy.  The  aimpoint  of  the  first  main  round  is  determined  by  supplying 
a  random  normal  deviate  representing  the  system  delivery  error  contribu¬ 
tion  to  the  aimpoint  probability  distribution  in  accordance  with  the 
following  equations.  This  is  illustrated  in  Figure  2-1. 


(2-1) 


(2-2) 


where 

XCl  =  horizontal  center  of  impact  for  main  rounds  (inches) 

R1  =  normal  random  deviate  (nondimensional ) 

ctxb  =  standard  deviation  of  variable  bias  in  horizontal 
direction  (inches) 

XC  =  intended  horizontal  aim  center  of  target  (inches) 

XB  =  horizontal  fixed  bias  for  main  rounds  (inches) 

YCl  =  vertical  center  of  impact  for  main  rounds  (inches) 


2-1 


(XC,  YC)  =  INTENDED  AIMPOINT 
(XC1,  YCl)  =  ACTUAL  AIMPOINT 
XB,  YB  =  FIXED  BIAS 

dxB,  dYg  -  VARIABLE  BIAS  SIGMA  TIMES  NORMAL  RANDOM  DEVIATE 

Figure  2-1.  Main  Round  System  Delivery  Errors 


R2  =  normal  random  deviate  (nondimensional ) 

OyB  “  vertical  standard  deviation  of  variable  bias  (inches) 

YC  =  intended  vertical  aim  center  of  target  (inches) 

YB  =  vertical  fixed  bias  for  main  rounds  (inches) 

The  impact  point  of  each  round  is  a  function  of  the  ballistic 
dispersion  which  is  defined  by  a  probability  distribution  about  the  center 
of  impact.  The  impact  point  of  this  round  is  determined  by  supplying  a 
random  normal  deviate  representing  the  random  delivery  error  contribution 
to  the  ballistic  dispersion  probability  distribution  in  accordance  with 
the  following  equations.  This  is  illustrated  in  Figure  2-2. 


(2-3) 
(2-4) 

where 

XI  =  horizontal  impact  point  of  main  round  (inches) 

R3  =  random  normal  deviate  (nondimensional) 
oxR  =  horizontal  standard  deviation  of  random  error  (inches) 

Oyi  =  horizontal  standard  deviation  of  lay  error  (inches) 

aL 

R4  =  random  normal  deviate  (nondimensional) 

Y1  =  vertical  impact  point  of  main  round  (inches) 

R5  =  random  normal  deviate  (nondimensional) 

OyB  =  vertical  standard  deviation  of  random  error  (inches) 

Oyi  =  vertical  standard  deviation  of  lay  error  (inches) 

R6  =  random  normal  deviate  (nondimensional) 

All  other  terms  have  been  previously  defined. 


XI  =  XCl  +  (R3)(aj(B)  + 


Y1  =  YCl  +  (R5)(ayB)  +  (ayL)(R6) 
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(XC,  YC)  =  INTENDED  AIMPOINT 
(XCl,  YCl)  =  ACTUAL  AIMPOINT 
(XI,  Yl)  =  IMPACT  POINT 

^XR  ^YR  '  ERRORS  SIGMA  TIMES  NORMAL  RANDOM  DEVIATE 

dxL  d^L  “  LAY  ERRORS  SIGMA  TIMES  NORMAL  RANDOM  DEVIATE 

Figure  2-2.  Main  Round  Random  Delivery  Errors 
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A  test  is  performed  to  ascertain  the  location  of  the  impact  point 
relative  to  the  target.  If  the  round  failed  to  hit  the  target  and  the 
miss  is  not  sensed,  or  if  the  resulting  impact  on  the  target  failed  to 
meet  requisite  kill  criteria,  another  round  is  fired,  the  impact  point  of 
which  is  negotiated  as  in  Equations  2-3  and  2-4.  If  the  miss  of  the 
round  is  sensed,  as  determined  by  drawing  a  suitable  uniform  random  number, 
the  aimpoint  of  the  subsequent  round  is  modified  by  supplying  a  random 
normal  deviate  to  the  following: 


(2-5) 


(2-6) 


where 


XC2  =  modified  horizontal  center  of  impact  (inches) 

R7  =  random  normal  deviate  (nondiraensional ) 

°XS  “  horizontal  standard  deviation  of  sensing  errors  (inches) 

YC2  =  modified  vertical  center  of  impact  (inches) 

R8  =  random  normal  deviate  (nondimensional ) 
aY5  =  vertical  standard  deviation  of  sensing  errors  (inches) 

All  other  terms  have  been  defined.  This  transform  is  illustrated  in  Figure  2-3 

Round  impact  simulation  continues  as  described  in  Equations  2-3  and 
2-4  except  that  the  lay  factor  is  eliminated  since  the  aimpoint  has  been 
modified  by  the  incorporation  of  the  effect  of  the  sensing  factor  and  that 
the  modified  aimpoint  values  replace  the  superseded  aimpoint  values. 

The  above  sequence  is  repeated  for  each  sample  in  determining  main 
round  target  hits.  Only  the  first  hit  is  recorded  for  each  sample  and 
it  updates  the  cumulative  target  hit  probabilities  of  only  those  round 
numbers  equal  to  or  greater  than  the  hitting  round  number.  After  the 
Monte  Carlo  sampling  exercise  has  been  completed,  these  cumulative  target 
hit  probabilities  are  divided  by  9,604  (the  total  number  of  Monte  Carlo 
trials)  yielding  the  mean  probability  of  hitting  the  target  for  a 
designated  number  of  rounds  fired. 
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(XI,  Yl)  =  PREVIOUS  IMPACT  POINT 

(XCl,  YCl)  =  PREVIOUS  AIMPOINT 

(XC2,  YC2)  =  MODIFIED  AIMPOINT 

dxs  dY5  »  SENSING  ERROR  SIGMA  TIMES  NORMAL  RANDOM  DEVIATE 

(X2,  Y2)  =  NEW  IMPACT  POINT 

^XR,  ^YR  "  'RANDOM  ERRORS  SIGMA  TIMES  NORMAL  RANDOM  DEVIATE 

Figure  2-3.  Delivery  Errors,  Sensed  Miss  Condition 


Spotting  Rifle  Round 


A  small  caliber  rifle  (e.g.,  0.50  caliber)  mounted  on  the  carriage 
of,  and  ballistically  matched  to,  a  major  artillery  piece  can  be  utilized 
as  a  round  delivery  assist  device  before  firing  a  main  round  projectile. 

The  procedure  previously  outlined  for  establishing  the  aim  and  impact 
points  of  the  main  round  is  repeated  for  the  spotting  rifle  rounds.  Once 
the  spotting  rifle  round  has  achieved  a  target  hit,  the  center  of  impact 
for  the  spotting  round  is  adjusted  by  subtracting  the  target  aimpoint 
coordinates.  Then  the  aimpoint  for  the  first  main  round  fired  is  computed 
in  accordance  with  the  following  equations,  and  is  illustrated  in  Figure  2-4. 


(2-7) 

(2-8) 


where 

XC3  =  adjusted  spotting  round  horizontal  center  of  impact  (inches) 

R9  =  normal  random  deviate  (nondimensional ) 

YC3  =  adjusted  spotting  round  vertical  center  of  impact  (inches) 

RIO  =  normal  random  deviate  (nondimensional) 

All  other  terms  have  been  previously  defined. 

The  computational  sequence  for  finding  the  main  round  impact  points 
continues  as  described  in  Equations  2-3  and  2-4. 
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(XC,  YC)  =  INTENDED  AIMPOINT 
(XCl,  YCl)  =  ACTUAL  AIMPOINT 

(XC3,  YC3)  =  ADJUSTED  SPOTTING  RIFLE  ROUND  AIMPOINT 
(XC4,  YC4)  =  FIRST  MAIN  ROUND  AIMPOINT 


igure  2-4.  Delivery  Errors,  Spotting  Round  Assist  Condition 


TIME-TO-KILL 


A  significant  parameter  which  is  monitored  throughout  the  simulation 
model  for  each  trial  is  the  time-to-kill.  In  addition  to  the  projectile 
time  of  flight,  the  time-to-kill  reflects  crew  proficiency  times.  The 
exponential  equation  quantifying  crew  proficiency  times  is  derived  using 
empirical  data  from  human  factors  studies  and  requires  Monte  Carlo  sampling 
for  its  solution.  All  kill  times  exceeding  120  seconds  are  ignored  in 
the  compilation  of  time-to-kill  data  for  each  Monte  Carlo  trial  since 
each  of  the  elements  of  the  time  Pk  array  represents  a  2-second  interval 
from  time  0  to  120  seconds. 

Firing  Time  Environmental  Factors 

The  estimation  of  the  elapsed  times  for  crews  to  fire  rounds  from 
their  main  armament  in  actual  combat  requires  more  than  simply  stating 
maximum  or  average  rates  of  fire  of  the  various  weapons  under  test  con¬ 
ditions  generally  tending  to  favor  rapid  fire.  The  objective  is  to 
estimate  for  each  weapon  of  interest  the  distribution  of  times,  not  only 
average  times,  that  crews  would  need  to  fire  the  weapon  in  a  firing 
engagement.  The  approach  representing  times  to  fire  aimed  rounds  in 
combat  is  based  on  the  following  basic  definitions  and  related 
assumpti ons . 

Firing  Engagement 

The  firing  engagement  is  defined  as  the  firing  of  one  or  more  rounds 
at  a  target.  Although  the  main  armament  may  be  fired  by  the  gunner  or 
the  commander  of  a  crew,  it  is  assumed  that  the  gunner  fires  the  weapon 
at  the  target  indicated  by  the  commander.  The  engagement  begins  when  the 
commander  issues  the  alert  element  "GUNNER"  of  the  initial  fire  command. 
Times  needed  by  the  commander  to  detect,  recognize,  and  identify  the 
target  are  excluded  from  the  engagement.  However,  any  time  that  is 
subsequently  needed  to  rotate  the  turret  or  to  find  the  target  is  part  of 
the  engagement. 

The  first  round  fired  at  any  target  is  normally  assumed  to  be  loaded 
before  the  engagement  begins  and  all  targets  fired  at  are  assumed  to  be 
seen  by  the  gunner.  Times  discussed  here  do  not  apply  to  combat  situa¬ 
tions  where  the  presence  of  a  target  is  suspected  but  its  actual  location 
is  not  known.  Subsequent  rounds  are  considered  fired  if  the  target  is 
not  killed  by  the  first  round.  Firing  of  rounds  may  stop  for  a  variety 
of  reasons,  such  as  the  target  being  killed  or  having  moved  to  a  new 
position  where  it  cannot  be  effectively  engaged. 
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Combat  Conditions 

Conditions  representative  of  a  combat  situation  in  which  a  tank,  for 
example,  would  fire  at  a  target  are  defined  as  combat  conditions.  Some  of 
the  factors  associated  with  combat  conditions  are: 

(1)  The  time  required  to  fire  the  first  round  at  a  target  depends 
on  whether  or  not  the  gunner  can  find  and  properly  aim  at  the 
target  indicated  by  the  commander.  Acquisition  of  the  target 
by  the  gunner  usually  presents  no  problem  under  demonstration 
or  test  situations,  but  in  the  field,  factors  such  as  obscura¬ 
tion  and  the  degree  of  contrast  between  the  target  and  its 
background  can  influence  the  time  required  by  the  gunner  to 
see  and  aim  at  the  target. 

(2)  The  gunner's  care  in  achieving  his  fine  lay  for  firing  is 
related  to  the  apparent  size  of  his  target.  The  degree  of 
care  that  he  chooses  to  exercise  in  making  a  fine  lay  for 
firing  especially  affects  the  first  round  firing  time. 

It  is  well  known  that  there  are  interactions  between  factors 
affecting  hit  probabilities  and  factors  influencing  firing  times.  Such 
interactions  obviously  must  be  taken  into  account  to  ensure  that  firing 
times  are  compatible  with  corresponding  weapon  delivery  accuracy  esti¬ 
mates.  From  the  viewpoint  of  how  quickly  a  tank  crew  can  fire  rounds, 
combat  conditions  are  categorized  in  the  following  three  ways: 

(1)  Favorable  conditions  are  combat  conditions  which  permit  all 
crews  to  perform  necessary  tasks  as  rapidly  as  they  can.  Time 
needed  by  various  crews  to  fire  rounds  still  varies  according 
to  their  proficiency.  Examples  of  factors  making  conditions 
favorable  would  be  the  absence  of  obscuration,  resulting  in 
the  gunner  experiencing  no  difficulty  in  finding  his  target, 
and  a  large  apparent  target  size,  making  the  exercise  of  great 
care  for  the  fine  lay  unnecessary. 

(2)  Intermediate  conditions  are  combat  conditions  midway  between 
favorable  and  severe.  Such  conditions  are  thought  of  as 
average. 

(3)  Severe  conditions  are  combat  conditions  which  limit  to  the 
maximum  extent  the  ability  of  all  crews  to  perform  necessary 
tasks.  Time  required  to  fire  rounds  still  varies  with  crew 
proficiency.  Obscuration  caused  by  smoke  and  dust  could  make 
the  gunner's  target  acquisition  task  very  difficult  and  a 
small  apparent  target  size  could  influence  the  gunner  to  make 
an  extremely  careful  fine  lay. 
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It  is  difficult  to  forecast  how  frequently  the  various  conditions 
would  occur  in  actual  combat.  Nevertheless,  favorable  conditions  would 
probably  occur  less  frequently  in  combat  than  in  the  tests  conducted  to 
obtain  firing  time  data  or  for  crew  qualification  purposes.  Unless 
particular  efforts  to  balance  combat  conditions  are  made,  these  tests 
tend  to  be  conducted  under  conditions  that  can  be  considered  favorable. 
However,  some  exercises  conducted  have  involved  conditions  that  could  be 
judged  intermediate  and  perhaps  border  on  severe.  Note  that,  since 
absolute  precision  is  not  involved  in  these  explanations,  informed  judg¬ 
ment  will  play  a  large  role  in  attempts  to  relate  particular  data  or 
performance  levels  to  an  associated  set  of  environmental  conditions. 

Crew  Proficiency 

Crew  proficiency  is  defined  as  the  level  of  proficiency  that  results 
from  the  tank  crew's  basic  ability  and  the  training  from  which  it  has 
benefited.  In  this  discussion,  proficiency  is  considered  exclusively 
from  the  viewpoint  of  its  influence  on  how  quickly  crews  can  fire  rounds. 
Crew  proficiency  is  categorized  in  the  following  ways: 

(1)  Best  crew  rating  denotes  the  maximum  proficiency  that  any  crew 
can  be  expected  to  attain.  Crews  in  this  category  are  the 
exceptionally  fast  crews. 

(2)  Good  crew  rating  denotes  a  level  of  proficiency  midway  between 
that  of  a  best  crew  and  an  acceptable  crew. 

(3)  Acceptable  crew  rating  denotes  the  minimum  level  of  proficiency 
that  all  crews  can  be  expected  to  attain  and  surpass. 

Here  again,  absolute  precision  is  not  possible  and  determinations 
associated  with  crew  proficiency  require  judgment.  The  proficiency 
associated  with  best  crews,  for  instance,  under  any  particular  environ¬ 
mental  conditions,  is  intended  to  correspond  to  a  level  achievable  only 
by  crews  with  maximum  inherent  ability  and  maximum  feasible  peacetime 
training.  Most  crews  cannot  achieve  this  status  no  matter  ;how  much 
training  they  might  undergo,  and  crews  with  the  greatest  ability  may  still 
fail  to  qualify  as  best  if  training  opportunities  sufficient  to  develop 
and  maintain  their  proficiency  are  not  provided.  The  level  of  performance 
associated  with  acceptable  crews,  again  under  particular  environmental 
conditions,  is  intended  to  be  a  lower  bound  to  acceptable  performance. 

This  status  is  achievable  by  crews  with  the  least  inherent  ability  if 
they  are  provided  adequate  training  opportunities.  Note  that,  since 
the  level  of  proficiency  demonstrated  by  any  crew  in  particular  environ¬ 
mental  circumstances  depends  on  both  ability  and  training,  either  factor 
can  compensate  for  failings. 
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Aimed  Round  Time  Factors 


The  program  also  utilizes  time  factors  for  aimed  rounds  including 
first  round  firing  time,  flight  time,  fixed  and  variable  times,  median 
times,  and  minimum  times  to  compile  time-to-kill  data.  These  factors 
are  discussed  in  the  paragraphs  that  follow. 

First  Round  Firing  Time 

This  is  defined  as  the  time  elapsing  between  target  recognition 
and  firing  of  the  first  round  assuming  that  the  decision  to  engage  the 
target  follows  target  recognition  without  any  delay.  Operations  involved 
in  firing  of  the  first  round  include,  as  applicable,  the  following: 

(1)  Target  recognition  by  the  gunner. 

(2)  Slewing  of  the  weapon. 

(3)  Ranging,  which  may  be  done  with  equipment  such  as  a  rangefinder, 
or  by  visual  estimation  of  the  range  to  the  target. 

(4)  Laying  of  the  weapon. 

Note  that  some  of  these  operations  can  be  carried  out  simultaneously. 

Loading  of  the  first  round  may  be  included  or  excluded.  In  the 
latter  case,  the  appropriate  type  of  round  is  considered  already  loaded 
before  recognition  of  the  target. 

Flight  Time 

Flight  time  is  the  time  that  elapses  between  the  time  a  round  begins 
to  be  propelled  forward  in  the  gun  and  the  time  it  reaches  the  range  of 
the  target. 

Subsequent  Round  Firing  Time 

This  is  defined  as  the  time  elapsing  between  firing  of  any  round 
and  firing  of  the  next  round  against  the  same  target.  Subsequent  round 
firing  time  involves,  as  applicable,  the  following: 

(1)  Flight  time  of  the  previous  round. 

(2)  Round  loading. 

(3)  Relaying  of  the  weapon,  which  may  or  may  not  reflect  efforts 
to  improve  the  aimpoint  with  reference  to  the  aimpoint  used 
for  the  previous  round. 
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Fixed  Times 

Fixed  time  is  defined  as  that  portion  of  total  firing  time  that 
is  independent  of  environmental  conditions,  crew  skill,  level  of  training, 
or  other  causes  of  variation  in  firing  time  performance.  This  time  is 
frequently  associated  with  mechanical  operations  such  as  automatic 
ramming.  Fixed  time  may  equal  zero. 

Variable  Times 

Variable  time  is  defined  as  that  component  of  total  first  or 
subsequent  round  firing  time  that  needs  to  be  represented  by  a  distribu¬ 
tion  of  times,  rather  than  by  a  fixed  time  component.  On  any  particular 
firing  occasion,  total  firing  time  equals  the  sum  of  the  applicable 
fixed  time  component  and  a  particular  time  from  the  distribution 
describing  the  variable  time  component.  Variable  time  is  tied  to  human 
operations  and  generally  will  not  equal  zero,  because  there  are  so  many 
causes  of  variation  that  are  not  controllable. 

Median  Times 

Median  times  of  the  distributions  apply  to  the  variable  component 
of  both  the  first  round  firing  time  and  the  subsequent  round  firing 
time.  Median  time  can  be  defined  as  the  particular  time  which  is  greater 
than  half  the  time  represented  by  the  corresponding  distribution  and 
smaller  than  the  remaining  half. 

For  the  special  case  where  a  logarithmic-normal  distribution  is 
fitted  to  data  consisting  of  N  particular  times,  the  median  time  is  the 
Wth  root  of  the  product  of  all  times  in  the  set  of  data  considered. 
Alternatively,  for  this  special  case,  the  median  time  corresponds  to 
the  logarithm  equaling  the  arithmetic  mean  of  the  N  logarithms  corre¬ 
sponding  to  the  times  constituting  the  basic  data.  Detailed  studies 
made  many  years  ago  in  the  U.S.  established  that  logarithmic-normal 
distributions  seemed  to  correspond  more  closely  than  other  known  distri¬ 
butions  to  test  data  obtained  for  tank  weapons  of  interest  at  the  time. 
Distributions  of  this  type  are  still  considered  useful  to  represent 
variable  components  of  firing  times  for  weapons  of  current  concern. 

Variability  Factors 

Variability  factors  indicate  to  what  extent  the  times  represented 
by  logarithmic-normal  distributions  vary  with  reference  to  the  applicable 
median  times.  A  standard  deviation  (based  on  calculations  with  natural 
logarithms)  of  about  0.5  has  been  found  to  apply  consistently  and  is 
normally  used,  unless  there  is  a  specific  reason  for  doing  otherwise. 

The  mathematical  discussion  that  follows  should  facilitate  under¬ 
standing  loqarithmic-normal  distributions  based  on  a  standard  deviation 
of  0.5. 
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Let  time  data  corresponding  to  a  particular  operation  and 

associated  conditions  be  v1 ,  v2,...,  vN  in  seconds.  Let  VI,  V2 . 

VN  be  logarithms  of  vl,  v2,...,  vN.  Assume  VI,  V2,...,  VN  are  normally 
distributed  with  median,  M,  and  standard  deviation,  S.  Let  m  =  antilog(M). 
Let  g(v)  =  cumulative  probability  that  a  particular  crew  performs  opera¬ 
tion  in  V  or  less  seconds.  Let  V  =  log  (v). 


/ML  1  19 

S  - - -  exp(- u^)  du  (2-9) 


where 

V  -  M  ^  log(v)  -  log(m)  ^  ioq[(l)l/^] 

s  s  ■' 


Everything  that  precedes  is  valid  for  both  natural  logarithms  and 
logarithms  to  the  base  10,  provided  one  or  the  other  is  consistently 
employed.  From  this  point  on,  only  natural  logarithms  will  be  considered 
because  computer  programs  developed  for  effectiveness  calculations  are 
based  on  their  use. 

The  particular  value  of  0.5  for  the  standard  deviation  S  is  used 
in  what  follows.  Consider  time  v  =  2m: 

=  Jin(4)  =  1.386 

g{v)  =  0.92 


Also,  it  is  easily  shown  that  if  v  =  m/2,  g(v)  =  0.08.  Consequently,  a 
general  rule  is  that  for  S  =  0.5,  all  but  about  1/6  of  the  times 
represented  by  a  logarithmic-normal  distribution  are  bracketed  between 
half  the  median  and  twice  the  median.  Likewise,  consider  time  v  =  3m: 


in  =  £n(9)  =  2.197 

g(v)  =  0.98 


where 

V  =  m/3,  g(v)  =  0.02 
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It  is  then  evident  that  all  but  about  1/25  of  the  times  in  a 
logarithmic-normal  distribution  for  which  S  =  0.5  lie  between  a  third 
of  the  median  and  three  times  the  median.  Similar  calculations  can  be 
made  to  help  one  understand  any  numerical  value  of  the  variability 
factor  S. 

Minimum  Times 

Minimum  time  is  defined  as  the  least  time  considered  possible.  Such 
a  least  time,  if  one  other  than  zero  is  specified,  overrides  any 
unrealistically  shorter  times  that  may  be  implied  by  a  particular  time 
distribution  used. 

The  data  base  used  to  establish  firing  time  estimates  for  various 
weapons,  yields  two  points  of  interest: 

(1)  Reliance  is  not  placed  on  any  one  source  of  data.  A  balanced 
viewpoint  stripped  of  biases  that  an  analyst  could  form  as 

a  result  of  specific  limited  experiences  is  required. 

(2)  Data  from  older  sources  need  not  be  discarded  as  newer  weapons 
become  of  principal  interest.  Several  relationships  can  be 
useful  even  after  weapons  used  to  obtain  underlying  data 
become  obsolete. 

As  in  the  case  of  delivery  accuracy  estimates,  the  construction  of 
firing  times  for  a  specific  weapon  requires  a  thorough  understanding 
and  definition  of  the  system  being  examined.  This  means  knowledge  of 
the  weapon  in  isolation  and  in  a  vehicle.  Times  of  flight,  slew  rates, 
target  designation  facilities,  fire  control  equipment,  standard  opera¬ 
ting  procedures,  breech  operation,  loading  devices,  round  weights,  and 
round  sizes  are  some  of  the  factors  which  can  affect  firing  time  esti¬ 
mates.  Tank  and  target  movement,  in  conjunction  with  the  assumed  condi¬ 
tions  and  system  definition,  also  affect  firing  times.  The  bulk  of 
information  available  on  firing  times  is  for  stationary  tanks  firing  at 
stationary  targets,  but  firing  time  estimates  should  be  made  for  other 
tank  target  movement  conditions  as  part  of  this  assessment. 

Main  Round 


At  the  outset  of  each  sample  simulation  run,  the  initialized  time 
of  zero  is  augmented  by  the  time  of  flight  of  the  main  round.  Henceforth, 
the  time  is  updated  by  crew  proficiency  factors  which  indicate  how  quickly 
crews  can  fire  rounds  as  a  result  of  their  training  and  basic  ability. 
Inherent  in  the  concept  of  crew  proficiency  would  be  such  elements  as 
target  acquisition  time,  aiming  time,  and  firing  time. 
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The  following  exponential  equation  derived  using  empirical  data  from 
human  factors  studies  requires  Monte  Carlo  sampling  for  the  solution 
of  a  time  factor  for  variable  crew  proficiency. 


HWM  =  (XM)  exp[(R9){ST)] 


(2-10) 


where 


HWM  =  incremental  time  factor  for  variable  crew  proficiency 
(seconds) 

XM  =  variable  median  time  (seconds) 

R9  =  random  normal  deviate  (nondimensional ) 

ST  =  slope  of  time  distribution 

If  HWM  is  less  than  a  minimum  reasonable  variable  time  to  accomplish  the 
firing,  the  predetermined  prescribed  minimum  time  to  fire  is  used  for 
HWM.  The  time  of  round  impact  is  then  computed  in  accordance  with  the 
equation: 


Tj  =  T^  +  HWM  +  TFS 


(2-11) 


where 


Tj  =  time  of  round  impact  (seconds) 

T^  =  current  accumulated  time  (seconds) 

TFS  =  fixed  median  time  (seconds) 

All  other  terms  have  been  defined. 

For  firing  subsequent  rounds,  the  sequence  returns  to  execute 
Equation  2-10  and  continues  as  described.  When  a  target  kill  has  been 
accomplished  for  a  given  kill  criterion,  the  time  of  kill  is  matched 
to  the  corresponding  time  interval  element  in  the  kill  array.  Only  the 
first  kill  is  recorded  for  each  sample  for  a  given  kill  criterion.  It 
updates  the  cumulative  target  kill  probabilities  of  only  those  times 
equal  to  or  greater  than  the  computed  kill  time.  After  the  Monte  Carlo 
sampling  exercise  has  been  completed,  each  cumulative  target  kill 
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probability  is  divided  by  9,604  (the  total  number  of  trials),  yielding  a 
mean  probability  of  killing  the  target  at  a  specific  time  interval  less 
than  or  equal  to  120  seconds. 

Spotting  Rifle  Round 

The  procedure  previously  outlined  for  establishing  the  time-to- 
kill  target  for  the  main  rounds  is  followed  for  the  spotting  rounds  to 
register  the  target  hit  time  which  is  then  added  to  the  flight  time  of 
the  first  main  round.  The  computational  sequence  for  finding  the  time 
to  target  kill  for  the  main  rounds  then  continues  as  described  in 
Equations  2-10  and  2-11. 

MEAN  VALUES 

The  mean  target  Pk  data,  the  mean  time  Pk  data,  the  mean  target  Ph 
data,  and  the  mean  time  Ph  data  are  computed  by  dividing  the  number  of 
samples  into  the  accumulated  element  values  of  the  applicable  arrays. 

The  arithmetic  mean  is  defined  as  the  sum  of  all  the  observations 
divided  by  the  total  number  of  observations. 


(2-12) 


where 


i  =  mean  value  of  n  observations 
£•  =  ith  observation 

n  =  total  number  of  observations 
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CASE  ROUNDS  FIRED 

The  average  number  of  rounds  fired  to  achieve  a  kill  for  each  case 
is  computed  in  accordance  with  the  equation: 


r 

AVI  =X)  i  (Pi  -  Pi.]) 
i=l 


(2-13) 


where 

AVI  =  average  number  of  rounds  fired  to  achieve  a  kill  for 
a  designated  case  (excluding  non-killing  trials) 

r  =  number  of  main  rounds  fired  per  sample 
=  mean  target  Pk  after  firing  i  rounds 

i  =  intermediate  factor  representing  the  number  of 
rounds  fired 

The  above  equation  is  applicable  only  if  a  kill  is  achieved  for  a 
designated  case  on  each  of  the  9,604  Monte  Carlo  trials.  Otherwise, 
computation  of  the  average  number  of  rounds  fired  to  achieve  a  kill 
for  a  given  case  is  accomplished  utilizing  not  only  those  rounds  fired 
for  non-killing  trials  but  also  an  additional  error  estimate  based  upon 
the  process  of  extrapolation  to  the  limit. 

In  addition  to  the  specified  number  of  rounds  fired  per  trial  as 
dictated  by  the  input  data,  the  effectiveness  of  five  succeeding  rounds 
is  estimated  based  upon  the  mean  results  of  the  last  five  rounds  fired 
for  all  the  sample  trials.  The  estimated  scaled  data  are  used  in  the 
computation  of  the  average  number  of  rounds  fired  per  case  in  accordance 
with  the  equation: 

^r+5 

AV2  =>]  i  (P^  -  P^.-,)  +  AVI 
i*r+l 

where 

AV2  =  rP^5  -  rP^  +  5P^5  -  (P,  t  P^^,  +  P^j  +  P^3  ♦  P^4)  +  AVI 
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/  1  ^  \ 

Introducing  the  scaling  factor  of  („ - p - I  =  SF  such  that 

V'^r  ■  / 

1  +  Pi.5(SF) 

P.  = -  for  r  <  i  <  r  +  5,  converting  terms,  and  letting 

^  1  +  SF 

Pp+5  =  1*0,  then: 

r-1 

5  +  (SF)  2  P^ 

AV2  =  r{l  -  P^)  +  5  -  - -  +  AVI 

''  (1  +  SF) 


AV2  =  r(l  -  P^)  + 


/ 

5(SF)  -  (SF)I  Pi 
'i=r-5 


(1  +  SF) 


+  AVI 


The  value  of  the  denominator  (1  +  SF)  is  assumed  to  be  unity,  therefore. 


(2-14) 


where 

AV2  =  average  number  of  rounds  fired  to  achieve  a  kill  for  a 
designated  case  (including  non-killing  trials) 

All  other  terms  have  been  previously  defined. 
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CLOSED  FORM  SOLUTION 

Although  the  basic  program  methodology  of  resolving  the  system  and 
random  errors  in  determining  projectile  impact  locations  is  accomplished 
by  Monte  Carlo  sampling,  a  closed  form  solution  of  a  single  main  round 
is  performed  for  each  new  impact  data  set  so  that  comparable  Monte  Carlo 
results  can  be  reviewed  for  agreement. 

In  describing  a  weapon's  accuracy  against  a  point  target,  or 
its  ability  to  hit  where  it  is  aimed,  it  is  generally  assumed  that  the 
errors  in  elevation  and  azimuth  are  distributed  according  to  the  normal 
or  Gaussian  law.  The  dispersion  of  a  fired  round  is  a  function  of  the 
ballistic  characteristics  of  the  particular  projectile  configuration 
and  of  the  meteorological  conditions  existing  at  the  time  and  place  of 
the  weapon  firing.  Dispersion  is  considered  a  random  error.  The  effect 
of  the  aiming  error  is  inherently  biased  in  some  direction  and  is 
regarded  as  a  systematic  error.  The  combined  distribution  representing 
the  distribution  of  impact  points  in  an  area  is  a  normal  distribution  in 
two  dimensions. 

The  assumption  that  the  distribution  of  impact  points  is  normal  is 
predicated  upon  the  reasonable  agreement  found  between  prediction  based 
upon  the  assumption  of  normality  and  the  data  of  actual  impact  distri¬ 
butions.  The  presumption  that  the  impact  distribution  is  generated  as 
the  sum  of  a  number  of  independent  random  variables  suggests  the 
applicability  of  the  central  limit  theorem  which  would  imply  that  the 
distribution  of  the  sum  would  be  approximately  normal. 

The  probability  of  a  single  main  round  hitting  the  target  and  the 
probability  of  a  single  main  round  hitting  the  target  given  a  sensed 
miss  are  computed  at  the  outset  of  the  JMEM/SS  Model  For  Direct  Fire 
Computer  Program  utilizing  the  error  function  associated  with  the  normal 
curve.  Since  the  weapon  errors  in  the  horizontal  and  vertical  directions 
are  distributed  normally  and  occur  independently,  the  resulting  distri¬ 
bution  of  errors  in  the  two  planes  is  a  bivariate  normal  distribution 
which  can  be  described  by  the  function: 


(2-15) 
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where 


PA^  =  probability  impact  point  will  lie  in  area 

oy  =  horizontal  standard  deviation  of  impact  error  for  a  single 
main  round  (inches) 

oy  =  vertical  standard  deviation  of  impact  error  for  a  single 
main  round  (inches) 

A^  =  ith  target  grid  cell  area 

X  =  represents  horizontal  (X)  target  coordinate 

Y  =  represents  vertical  (Y)  target  coordinate 

Equation  2-15  is  used  to  compute  the  target  cell  hit  probabilities  for  a 
single  main  round  given  a  sensed  miss.  The  computation  can  be  used  for 
either  a  single  main  round  or,  if  the  option  is  exercised,  for  a  single 
main  round  with  a  spotting  rifle  round  assist.  The  hit  probabilities  for 
all  the  target  grid  cells  are  summed  to  include  the  entire  target  form  in 
accordance  with  the  equation: 


(2-16) 


where 

Ph  =  probability  of  hitting  target  with  a  single  round 
a  =  number  of  target  cells 
All  other  terms  have  been  previously  defined. 

In  order  to  compute  the  sum  of  standard  deviations,  the  root-sum- 
square  combination  is  executed  to  determine  the  standard  deviation  of 
impacting  errors  in  the  horizontal  and  vertical  directions  for  a  single 
main  round,  for  a  single  main  round  given  a  sensed  miss,  and  for  a  single 
main  round  with  a  spotting  round  assist. 
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Single  Main  Round 

Assuming  that  the  random  errors,  variable  bias,  and  lay  errors  for 
a  single  main  round  are  independent  and  normally  distributed. 


II 

X 

ID 

(2.2.2 
'’XB  °XL 

II 

>- 

D 

(2.2.2 
®YR  '^YB  '^YL 
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where 

Own  =  standard  deviation  of  horizontal  random  error  for  main 
rounds  (inches) 

CTvg  =  standard  deviation  of  horizontal  variable  bias  for  main 
^  rounds  (inches) 

ow,  =  standard  deviation  of  horizontal  lay  error  for 
main  rounds  (inches) 

Own  =  standard  deviation  of  vertical  random  error  for  main 
rounds  (inches) 

oyD  =  standard  deviation  of  vertical  variable  bias  for  main 
rounds  (inches) 

ayi  =  standard  deviation  of  vertical  lay  error  for  main 
rounds  (inches) 

All  other  terms  have  been  previously  defined. 
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Single  Main  Round  Given  Sensed  Miss 

Assuming  that  the  sensing  errors  and  random  errors  for  a  single 
main  round  are  independent  and  normally  distributed: 


°A  =  V 

^  2  ^  x2 

°XS  ^^®XR^ 

‘'B  =V 

^  2  .  ./  ^2 
ays  +  2(aYp) 

(2-19) 
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where 


an  =  horizontal  standard  deviation  of  impact  error  for  a  single 
main  round  given  a  sensed  miss  (inches) 

®XS  ~  standard  deviation  of  horizontal  sensing  error  for  main 
rounds  (inches) 

an  =  vertical  standard  deviation  of  impact  error  for  a  single 
main  round  given  a  sensed  miss  (inches) 

aye  =  standard  deviation  of  vertical  sensing  error  for  main 
rounds  (inches) 

All  other  terms  have  been  defined. 

Single  Main  Round  with  a  Spotting  Round  Assist 

Assuming  main  round  impact  sigma  and  the  sensing  errors,  lay  errors, 
and  variable  bias  from  a  spotting  rifle  round  are  independent  and  normally 
distributed: 


axi  "  +  [.5  (aj^5^  +  OxLl^^  ° 


2 

XRl 


Y1  V°Y  ^  ^®YS1  *  *^¥11^^^  '^YRl 
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where 


Oy,  =  horizontal  standard  deviation  of  impact  error  for  a  single 
main  round  with  spotting  round  assist  (inches) 

=  standard  deviation  of  horizontal  sensing  error  for  spotting 
^  rounds  (inches) 

o  =  horizontal  standard  deviation  of  lay  error  for  spotting 
rounds  (inches) 

ayp-,  =  standard  deviation  of  horizontal  random  error  for  spotting 
rounds  (inches) 

=  vertical  standard  deviation  of  impact  error  (inches) 

„  =  vertical  standard  deviation  of  vertical  sensing  error  for 

spotting  rounds  (inches) 

ay,  =  vertical  standard  deviation  of  lay  error  for  spotting 
rounds  (inches) 

oyon  =  standard  deviation  of  vertical  random  error  for  spotting 
rounds  (inches) 

All  other  terms  have  been  previously  defined. 
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SECTION  III 


CONCEPTUAL  FLOWCHARTS 


This  section  contains  conceptual  flowcharts  for  the  JMEM/SS  Model 
for  Direct  Fire  Computer  Program.  A  flowchart  for  each  significant 
section  follows  a  narrative  keyed  to  the  flowchart  by  step  numbers. 

Steps  are  indicated  on  the  flowcharts  by  numbers  enclosed  in  hexagons. 

These  flowcharts  provide  an  effective  pictorial  summary  of  the 
computational  steps  of  this  program.  A  more  detailed  level  of  discussion 
is  provided  in  Section  IV,  Simulation  Model. 


MAIN  ROUTINE 

The  MAIN  Routine  is  used  to  locate  the  impact  points  of  the  main 
and  spotting  rifle  rounds  and  to  compute  the  target  probability  of  kill, 
the  target  probability  of  hit,  the  time  probability  of  kill,  and  the 
time  probability  of  hit  data. 

The  following  description  is  an  outline  of  the  flow  of  the  MAIN 
Routine.  The  description  is  keyed  to  the  conceptual  flowchart,  Figure 
3-1,  by  the  numbers  enclosed  in  hexagons. 


(  1  /  Execution  is  begun  by  calling  Subroutine  INPUT  to  read  and 
write  the  input  data  for  the  9,604  Monte  Carlo  trials.  Sub¬ 
routine  SSHOT  is  called  to  compute  and  write  a  closed  form 
solution  of  a  single  main  round  hit  probability  and  of  a 
single  main  round  hit  probability  given  a  sensed  miss  before 
initializing  selected  arrays  in  preparation  for  later 
processing. 

/z)  After  the  Monte  Carlo  trial  loop  is  entered.  Subroutine 
' — '  NRAN31  is  called  to  provide  the  random  normal  deviates 

necessary  to  compute  the  impact  points  and  to  update  the 
time-to-kill  for  either  the  spotting  round  or  the  main  round. 
If  the  target  is  not  hit,  the  next  round  is  fired,  or  the 
sample  loop  is  re-entered  if  there  are  any  remaining  Monte 
Carlo  trials  to  be  negotiated.  If  a  target  hit  has  been 
achieved  by  a  spotting  rifle  round,  the  main  round  loop  is 
entered. 


The  kill  criteria  loop  for  the  main  round  is  entered.  When 
the  target  hit  achieves  the  first  target  kill  of  each  sample 
for  the  criterion  under  consideration,  the  cumulative  target 
Pk,  and  the  cumulative  time  Pk  arrays  are  updated. 
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If  there  are  no  more  rounds  to  be  fired  for  a  sample  and  all 
the  Monte  Carlo  trials  have  been  completed,  the  mean  values 
of  the  cumulative  target  Pk  and  cumulative  time  Pk  arrays  are 
computed  and  written. 

If  the  number  of  main  rounds  fired  per  sample  is  not  equal  to 
one.  Subroutine  KILRAT  is  called  to  compute  and  write  the 
average  number  of  main  rounds  fired  per  case.  Upon  return 
from  Subroutine  KILRAT,  either  a  new  input  data  set  is  read 
or  program  execution  terminates. 


Figure  3-1,  Flowchart,  Main  Routine 
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SUBROUTINE  KILRAT 


This  subroutine  is  used  to  compute  and  write  the  average  number  of 
main  rounds  fired  for  each  case  (which  represents  each  given  kill 
criterion  and  the  target  hit  data).  The  following  description  is  an 
outline  of  the  flow  of  Subroutine  KILRAT.  The  description  is  keyed 
to  the  conceptual  flowchart.  Figure  3-2,  by  the  numbers  enclosed  in 
hexagons. 

The  case  loop  is  entered  and  the  round  loop  is  executed  to 
compute  the  average  number  of  rounds  fired  for  a  given  case. 


<D 


If  a  kill  occurred  within  the  set  number  of  rounds  fired  for 
a  given  case  for  each  sample,  the  average  number  of  rounds 
fired  for  that  case  is  written.  Otherwise,  the  computation 
of  the  average  number  of  rounds  fired  for  that  case,  includ¬ 
ing  not  only  those  rounds  fired  on  non-killing  trials  but 
also  an  estimated  extrapolated  factor,  is  performed  before 
writing  the  results.  Control  returns  to  the  calling  routine 
after  completion  of  all  the  cases. 
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ENTER 


Figure  3-2.  Flowchart,  Subroutine  KILRAT 


SUBROUTINE  SSHOT  (PCH,  PCHl ) 


This  subroutine  is  used  to  compute  the  probability  of  hitting  the 
target  for  a  single  main  round  and  the  probability  of  hitting  the 
target  for  a  single  main  round  given  a  sensed  miss  by  using  a  closed 
form  solution.  The  following  description  is  an  outline  of  the  flow  of 
Subroutine  SSHOT.  The  description  is  keyed  to  the  conceptual  flowchart. 
Figure  3-3,  by  the  numbers  enclosed  in  hexagons. 

{T)  The  standard  deviation  of  the  impact  error  for  a  main  round 
and  the  standard  deviation  of  impact  error  for  a  main  round 
given  a  sensed  miss  are  computed.  If  spotting  rifle  round 
data  is  to  be  included,  these  standard  deviations  are  appro¬ 
priately  modified. 


All  grid  cells  describing  the  target  form  are  identified  in 
order  to  compute  and  to  appropriately  sum  the  probabilities 
of  a  single  round  hitting  these  designated  target  cells  by 
using  the  cumulative  normal  distribution  function. 
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Figure  3-3.  Flowchart,  Subroutine  SSHOT 

Next  page  is  blank. 
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SECTION  IV 


SIMULATION  MODEL 


The  JMEM/SS  Model  for  Direct  Fire  Computer  Program  is  comprised  of 
a  MAIN  Routine  and  six  subprograms.  This  section  describes  how  each  of 
these  elements  of  the  program  is  used  to  implement  the  Mathematical 
Model  described  in  Section  II.  The  subprograms  are  described  one  at 
a  time;  the  purpose  of  each  subprogram  is  defined,  and  the  function  of 
each  group  of  statements  within  the  subprogram  is  described. 

After  the  initial  data  set  pertaining  to  all  of  the  Monte  Carlo 
sample  trials  is  read  by  Subroutine  INPUT,  the  target  grid  and  Pk  arrays 
are  set  up.  Subroutine  SSHOT  is  called  to  compute  a  closed  form  solution 
of  the  probability  of  hitting  the  target  for  a  single  main  round  and  the 
and  the  probability  of  hitting  the  target  for  a  single  main  round  given 
a  sensed  miss. 

Subroutine  NRAN31  and  Function  RANDOM  are  utilized  in  applying  a 
Monte  Carlo  sampling  technique  in  determining  target  kills,  round 
aimpoints,  round  impact  points,  and  the  time-to-kill  values.  The  results 
are  generated  in  the  MAIN  Routine  from  9,604  Monte  Carlo  trials  which 
yield  mean  target  Pk  data,  mean  time  Pk  data,  mean  target  Ph  data,  mean 
time  Ph  data,  and  the  average  number  of  rounds  fired  per  case. 

In  the  following  paragraphs,  the  individual  FORTRAN  statements  that 
make  up  the  Simulation  Model  include  references  to  Section  II,  Mathe¬ 
matical  Model.  These  references  consist  primarily  of  the  boxed  equations 
that  gave  rise  to  the  FORTRAN  statements,  but  not  their  derivations.  The 
same  boxed  equations  can  be  found  in  the  Mathematical  Model  with  their 
derivations.  The  number  enclosed  in  parentheses  next  to  each  equation 
is  its  identifying  number  in  Section  II.  In  addition,  comment  statements 
appearing  in  the  source  listing  are  included  in  the  Simulation  Model  for 
the  sake  of  completeness. 

Several  subprograms  call,  COMMON  and  DIMENSION  FORTRAN  statements, 
whose  purpose  is  to  set  aside  arrays  to  store  values  of  the  variables 
appearing  in  the  statements.  The  definition  of  the  variable  names,  as 
well  as  the  variables  in  the  subroutine  call  statements,  are  provided 
in  the  List  of  Abbreviations  and  Symbols  at  the  end  of  this  section. 


n  r>  r> 


MAIN  ROUTINE  (1  of  13) 

This  program  is  used  to  locate  the  impact  points  of  the  main  and 
spotting  rifle  rounds,  to  monitor  the  time-to-kill ,  and  to  compute^ 
target  probability  of  kill,  target  probability  of  hit,  time  probability 
of  kill,  and  time  probability  of  hit  data. 

The  program  is  begun  by  execution  of  the  statements 
c 

c  jMFM/ss  Mnnn.  for  dirfct  ftrf 
c 

C  main  ppnr,RAM 
C  PURpnsp 

C  tnrATFS  IMPACT  POTMT3  OF  MAIN  AND  SPOTTING  RTFLF  ROIINDSJ  MONITORS 

C  THF  TTMF  TO  KTI.LJ  COMPUTFS  PK  AMD  PH  TARGFT  DATA  AND  PK  AND  PH  TIMF  DATA 

^  *  *  *************** 
C 

c  shrprogpams  rfdutrfd 

C  INPUT-  RFADS  INPUT  DATA, 

C  KTIRAT-  rn^PUTFS  AND  WRITFS  THF  AVFRAGF  NU'^BFR  OF  MAIN  ROUNDS  FIRED 
C  PFR  CASF. 

C  NRAM^I-  PROVIDES  NORMAL  DEVIATES  FOR  MOMTE  CARLO  SAMPLING 
C  R  ANDnM-HNTFf)RM  RANDOM  WIJMRFR  GFNFRATOR. 

C  SSHOT-  COMPilTFS  PH  FO'?  A  SINGLE  main  ROUND  AND  PH  FDR  A  SINGLE  ROUND 
C  GIVEN  A  SENSED  MTSS  USING  A  CLOSED  FORM  SOLUTION, 

C 

C*****<k***********************************^*************^*********************** 

c 

dimension  PKI^OOO  ,P)#AVfP0tIiAY(S1  ) 

rOMMON  PK. AX, AY ,SIGXR,STGYR, STGXR,  S I G Y B , S T GX S • S T G Y S , PROPS , X B # YR , 

1  xr  .  YC  .REI.  .TFl.TFS.XMi,  X  Mp ,  y  OF  F  •  •'»OP  T  I  ,  NROS »  NX  ,  NY  ,  FL  T  •  S  T  1  ,  ST2,  AMJ  1  , 

PAmtp,  TK  (Sn.pnoi  ,  7. 71  ,  AK  TL  (^>1  .  RI  •  STGXI  ,  STGYL  ,RELF 

^.Sir>XS1.STGYS1.STGXBt,STGYB1,STGXRt,STGYR1,PRSt»XRl,YB!,Fl.  T1,TFll# 
alESI  ,  XMfi  , XMpi , AMTii . AMTPt .ST  n . ST?1  rSTGXLl #STGYL I #NRDS1 .PASSN 

common/  randm/ti 

c 

which  are  used  to  allocate  storage  for  arrays  and  to  provide  for  the 
exchange  of  data  between  the  MAIN  Routine  and  the  remainder  of  the 
program.  The  random  number  seed  is  initialized,  and  Subroutine  INPUT 
is  called  to  read  and  write  the  input  data  by  execution  of  the  statements 


c  SET  71=  TO  ANY  ODD  NUMBER  FOR  RANDOM  SEED, 
C 

11=11111111 
READ  INPUT 
?  CAl  L  INPUT 


The  statements 


C  rOMPilTF  PPn^AHTlITV  OF  HTT  (IF  FIRST  ROUND  AND  PRODAnil  TTV  OF  MIT  GIVFM 
c  A  SFNSFD  MISS  IN  A  CLOSFD  FORM  SOLUTION  AND  WRITF  RFSULTS. 

C 


4-2 


MAIN  ROUTINE  (2  of  13) 

r Ai  I  ssHfUCPCM.prMt  i 
PfH.Pr.HI 

hO  Fn«MATf/^1H  PPOH.  OF  HTT  np  FTP<5T  R^U^m  =  ,F16.7,/?ftH  PROR.  OF  HIT 
l/SFMSFD  WTSS  *  ,F1h.71 

are  used  to  call  Subroutine  SSHOT  in  order  to  compute  a  closed  form 
solution,  i.e.  without  utilizing  a  Monte  Carlo  sampling  technique,  of 
a  single  main  round  target  hit  probability  and  of  a  single  main  round 
target  hit  probability  given  a  sensed  miss.  These  results,  which  should 
be  comparable  to  subsequent  Monte  Carlo  derived  values  for  the  same 
target/weapon  parameters,  are  written  before  the  Monte  Carlo  sample  size, 
the  time-to-kill  array,  and  the  mean  Pk  array  are  initialized  by  the 
execution  of  the  statements 


c 

C  INITI4LT7F  SAMPlF  St7F. 

C 

c 

C  TNTTTAtT7F  T  I 'IF  ARRAY  AND  PPHRARTlITY  AWRAY. 
C 

nORS.Isl  0h 
nniSMTri 

1SR  AKTl  f T.JlsO. 

nn  QS  Tsi .NROS 
qs  7fT.J)sO. 


The  statements 


c 

c  PRTMT  nuT  SAMPI.f  ST7F. 

n 

WRTTF  f^.^PUlSAMP 

C 

C  start  SAMPI  TNf* 

c 

THl ?=0 
THI 

m  qq  Tsi . nsamp 
C 

C  TNTTTAI.T7F  VARTAflI.FS  FOR  FACH  SAMPt.F. 
C 

SHMCAS  sft. 
yri sn, 
vriso. 

THtOsO 
THI  t  =0 
T  =  n. 

TH^sO 
ooqujri 
qq  rr.nsn. 


are  used  to  write  the  sample  size,  to  enter  the  loop  of  the  9,604  Monte 
Carlo  samples,  and  to  initialize  the  following  variables  for  each  sample: 
the  personnel  casualty  counter,  the  center  of  impact  of  the  main  or 
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nor)  non  non 


MAIN  ROUTINE  (3  of  13) 


spotting  rounds,  the  main  round  flag,  the  spotting  rifle  round  counter, 
the  time-to-kill ,  the  first  main  round  hit  flag,  and  the  kill  flag  array. 

The  statements 

TF  nMLV  main  rounds  CONSIOFRfO  GO  TO  201. 

TF  fNnPTr.NF.?1GOT020I 
rilL  NRANSI (R1 .R?) 

SFT  T  TO  TTMF  f)F  FLIGHT  OF  ONF  SPOTTING  RIFLF  SHOT. 

TsFITI 

XC1.TC1  IS  fFNTFR  OF  TNPAtT  FOR  FIRST  ROUND  OF  SPOTTING  RIFIF  SHOTS'. 

yr  1  sRI  *STG)(H1  ♦XC+KRI 
vr 1 =R?*STGVB1 fYC+THI 


are  used  to  check  the  flag  governing  program  options.  If  NOPTl  does  not 
equal  two,  the  program  considers  only  main  rounds  and  control  is  trans¬ 
ferred  to  Statement  201.  If  NOPTl  equals  two,  the  time-to-kill  is  updated 
by  the  time  of  flight  of  a  spotting  rifle  round,  and  Subroutine  NRAN31 
is  called  to  provide  the  random  normal  deviates  needed  to  compute  the 
horizontal  and  vertical  coordinates  of  the  center  of  impact  for  the  first 
spotting  rifle  round. 

The  statements 


FACsI. 

C 

C  NOPMAI  OFVTATFS  FOP  RAMDOW  AMO  LAY  FRROPS  OF  SPOTTIMG  RTFLF. 

C 

200  CALL  NRAM-^1  f  PI  ,  W?) 

TALL  NRAN-^1  (R 
C 

C  THtt  TOFNTTFTFS  THF  SPOTTING  RTFIF  ROUNOS. 

C 

TH1 1=TH1 1 ♦I 
C 

C  ORAW  NORMAI  nFVTATFS  -OR  TT^^F  FOR  SPOTTING  RTFLF. 

C 

TALL  NPAN-^1  (T1  . 

r 

C  FOR  FIRST  ROIJNO  11*^?  OF  SPOTTING  RIFLF  GO  TO  ^10. 

C 

TFT  TH1  1  .FO.  1  IGOTflP!  0 
C 

C  FOP  SUflSFQMFNT  ROUND  TTMfS* 

C 

H 1  ♦P  YP  (  T  1  *STP  n 

TF  f  HWM.I  F,  A^TPI  )  Ha/M  =  AMT?1 

TsT+HaM+'TFSI  • 
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MAIN  ROUTINE  (4  of  13) 

RHTn^l 1 

?10  HWMsX^i t *FXPf T1 1 n 

TF  fHwM.I  F. AVT1  1  )  M^MrAMTl 1 
TsTfH.^M  +  TFII 
?11  rnMTTNMF 

c 

C  TMPATT  PnTNT  FOR  THIS  SPOTTING  RTFI.F  SHOT. 

C 

y  1  sXCI  *STGXR1  t*^^*??TGXI.  1  *F  AC 
YlsYriFR?ASTGVRl^Ril»STGVl  IaFAC 
C 

C  GO  SFF  TF  THIS  SHOT  HTT. 

C 

GOTOlhO 

are  used  to  compute  the  impact  point  and  to  update  the  time-to-kill  the 
target  for  the  spotting  rifle  rounds  after  calling  Subroutine  NRAN31  in 
order  to  furnish  the  random  normal  deviates  for  the  computations.  Con¬ 
trol  is  then  transferred  to  Statement  160  to  ascertain  if  the  spotting 
round  impact  point  lies  within  the  extremities  of  the  target  area. 

The  statements 


c 

C  DRA^  RANDOM  DFVTATFS  FOR  VARIARLF  BIAS  FRRORS  FOR  FTRST  MAIN  POUND. 
C 

?01  CALL  NRAN^I fPI .R?) 

C 

C  UPDATE  TTmf  RY  FLIGHT  TTMF  OF  MATN  ROUND. 

C 

T=TfFlT 

C 

C  UPDATE  CFNTFl^  OF  IMPACT  FO»  MATN  ROUNDS. 

C 

XC1 =R1 *STGXfl+XCFXR+XC1 
YCI rRP*STGYR+YC^YR*yCt 
C 

C  FAC  TS  A  FACTOR  TO  TNCL JDF  LAY  FRROR  TF  IT  IS  1. 

C 

FAC=1. 


are  used  to  set  the  lay  error  flag  and  to  update  the  time-to-kill  and 
center  of  impact  for  the  first  main  round  of  each  sample  after  calling 
Subroutine  NRAN31  in  order  to  provide  the  random  normal  deviates  required 
to  compute  the  horizontal  and  vertical  aim  coordinates  of  the  first  main 
round  in  accordance  with  the  following  equations. 


First  main  round: 


XCl  =  {R1)(oxb)  +  XC  +  XB 


(2-1) 
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YCl  =  (R2)(CTYg)  +  YC  +  YB 


First  main  round  subsequent  to  spotting  rifle  round: 


(2-2) 


(2-7) 


(2-8) 


FNTF«  matM  WOUlViD  Lonp. 

.1  =  0 

PAS  J  =  J*1 

THtn  TS  A  SWITCH  in  TFIJ.  IF  THIS  TS  MAIN  RDUMD  OP  *  TH10=1  TS  MAIM 
ROUf'JO. 

THIOrl 
TH?0  =  0 

SAMPt  F  IMPACT  POTNI  FOR  IHTS  ROUMF) 

IF  fF  AC  .FO.  0 . 1  Cnim  0  1 

IF  RDlJNn  HAS  MFRM  SFMSpn  PFFORF,  F^CliJOF  IMF  LAY  FRROR. 

CAM  ^RAM^I f R1 ,R?) 

XP=STGXL*Ri 

VP»STfiYL*R? 


are  used  to  initialize  the  main  rounds  fired  counter,  to  execute  the 
main  round  loop,  to  set  the  main  round  flag,  and  to  test  the  lay  error 
flag.  If  a  main  round  has  not  been  sensed  previously.  Subroutine  NRAN31 
is  called  to  supply  the  random  normal  deviates  required  to  quantify  the 
lay  error  in  the  horizontal  and  vertical  directions. 

The  statements 

tot  CALL  NRAN^i fRl .R?! 

C 

C  FtNO  IMPACT  PDTWT  OF  THIS  '^ATM  ROHNO. 
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MAIN  ROUTINE  (6  of  13) 
c 

xisxri4.Pi*STr,)(»^FAr*x? 

Y1  svn  ♦P?*STr,YRfFAC*YP 

IJPOATF  THF  TT'^F  FOP  THTS  '^ATN  PIU’NJD* 

TA!  I,  NRAN^I  f  T1  ,7?) 

TFf .f.FQ.  t  )anTr!i SR 
HWMsX^pAf XP f T 1 *ST?1 
TF  (H'Vm'.i  F,  HWVsAMT? 

T=TfHv«M*TFS 
nnrni ho 

1S9  HW^sxMi  *FyP(T1  *STn 

IFfNiNM.i  F.  AMTI)  HWMrAMTI 
T=T>MWM+TF 1 


are  used  to  update  the  time-to-kill  the  target  and  to  compute  the  main 
round  impact  point  after  calling  Subroutine  NRAN31  which  provides  the 
random  normal  deviates  necessary  for  these  computations  in  accordance 
with  the  equations: 


XI  =  XCl  +  (R3)(aj(^)  + 


Y1  =  YCl  +  (R5)(aY[^)  +  (aYL)(R6) 


HWM  =  (XM)  exp[(R9)(ST)] 


Tj  =  1^+  HWM  +  TFS 


(2-3) 

(2-4) 

(2-10) 

(2-11) 


If  the  target  is  defiladed  on  any  main  or  spotting  round,  or  if 
the  impact  point  of  any  main  or  spotting  round  lies  outside  of  the 
enveloping  target  grid  cells,  or  if  the  identifying  kill  array  has  no 
target  kill  criteria  associated  with  the  cell  which  is  identified  as 
hit,  control  is  transferred  to  Statement  1  by  execution  of  the  statements 


IS  TA»GFT  OFFILAGFO? 

IfcO  TFfYt.lT’.YOFF)  GO  ID  1 


C 

C 

C 


MAIN  ROUTINE  (7  of  13) 
c 

C  TS  THTS  PnitND  WTTHTM  VULNFR&R  T  L  T  T  Y  MATRIX? 

C  .  , 

TFf XI AX  r 1  ) IGHTOI 
TFf  XI  .r,T.  AXfMXf  n  IGOTOI 
TFfYi^i.T^AYf  1  )ir,nTni 
TFCYI  .r.T.AY(NYf1  DGOTOl 

C 

C  WHFRF  TS  THIS  ROlINO  TN  THF  VUL  NFR  AR  TL  T  T  Y  MATRIX? 
C 

.M  =  1 

5  Jlsllfl 

IF  f  Yi  ,r*T‘.AYr  Ji  )  inoTn-^ 

51  .!?=1 

<1  .T?z|?+1 

TFf  <1  ,nT.Ax(,i?i  innina 

c 

C  TS  THTS  ROliNO  A  HIT? 

C 

TF  r  TK  f  Jt  .  .T?T  .L  F*0  TROTni 


Then  the  statements 
c 

C  TF  THTS  TS  A  SPOTTING  RTFIF  HIT  M^OATF  MLIMRFR  OF  pnijNOS  TT  TOOK  AMO 
C  START  F  TRTMG  '  THF  ma  T  N  ROtlND. 

C 

TFfTM10/.Fr3.0)  THfi=THli+THt1 
TFf TMtO.FQ.ninOTOPOU 

are  used  to  update  the  total  number  of  spotting  rounds  fired  to  achieve 
a  spotting  rifle  round  target  hit  before  transferring  control  to  State- ^ 
ment  204  to  fire  a  main  round.  Diminution  of  weapon  effectiveness  attri¬ 
buted  to  the  reliability  of  the  round  is  achieved  by  the  execution  of  the 
statement 

c 

C  TS  THTS  HIT  A  RFLTARLF  ROmnD? 

C 

TFfRAMnOMf  t 1 .GT ,RFL TGOTOR? 


The  statements 


c 

C  TF  TH-TS  HIT  TS  NOT  THF  FIRST  MTT  GO  TO  SS 

c 

TFf THw.GT.niGnTOSS 
C 

C  UPRATF  THF  PRRA-RTLTTY  ARRAY  AND  T  T  MF  ARRAY  FOR  THTS  HTT* 
C 

THWsI 

.nos. I 

nnsiKs.no^NRDS 
51  7rK,Ms;fK.6)  +  i . 

,11  OsP.4T/P. 
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MAIN  ROUTINE  (8  of  13) 

TFf Ji o.HT.hi inninss 
OntSlJl 1sJ10,61 

1S1  AKTl  f Ji  1  .MaAKTK  T1  1  +  1  . 

are  used  to  update  the  cumulative  results  of  the  target  and  time  arrays 
by  identifying  the  first  main  round  hit  and  time  for  each  sample  case. 

The  statements 


C  nOFS  THTS  HTT  HAVF  A  PFLTArtl.F  FU7F? 

C 

SS  TFfRAMnnMd ).GT,9FLF)nnTOQ7 
TMPOst 

are  used  to  check  fuze  actuation  and  to  set  the  fuze  actuation  flag. 
The  statements 


ARF  r.nMSIDFRFO 


c 

C  TS  THIS  HTT  A  KTIL? 

r 

HWPrWANOnMf 1  ) 

KPP  =  TK(  M  ..l?T 
MW  1=0* 

c 

C  DTFFFRFNT  KILL  CRTTFPTA 

C 

nOP^KPtrl .d 

TFftfKPI T.GT.O.TGnTDQ^ 

TFT PK ) .LT.HWPlGnTOR? 

C 

C  TF  THTS  wtT  TS  A  KTLL  FOR  THTS 
C  TTMF  ARRAY. 

C 

,T10  =  P.^T/P. 
TFfJiO.GT.bt^GOTDIS? 
nnt  s^.n  i  =.ti  o.6i 

1S^  AKTI  f.II  1  .K?n  =  AKTI.  CJn  .K?n  +  1 
IS?  nna?K?n=  f .nros 

n?  7fK?0,K?n=7fK?0.K?tTfl. 
rf K?i >=i . 

GfTTOR^ 

9?  MWlzt. 

9^  rOMTMlIF 


CRTTFRIOM  HPOATF  PRORARILTIY  ARRAY  AND 


are  used  to  set  the  number  of  kill  criteria  and  to  execute  the  kill 
criteria  loop.  If  the  kill  flag  array  element  designating  a  given  kill 
criterion  is  greater  than  zero,  (indicating  achievement  of  target  kill 
for  that  criterion  by  a  previous  round)  the  loop  is  re-entered  for 
execution.  If  the  Pk  value  for  the  designated  kill  criterion  of  the 
impacted  target  grid  cell  is  less  than  a  selected  uniform  random  number 
between  zero  and  one,  the  flag  indicating  failure  to  achieve  all  kill 
criteria  is  set.  If  the  target  hit  achieves  a  kill  for  the  criterion 
under  consideration,  the  probability  of  target  kill  array  and  the  time- 
to-kill  array  (providing  that  the  time-to-kill  does  not  exceed  120 
seconds)  are  updated  to  ultimately  reflect  the  cumulative  results  of 
all  9,604  sample  cases. 
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MAIN  ROUTINE  (9  of  13) 


The  statements 
c 

C  TT  THIS  ROUND  “'AS  NOT  A  KTIK  FOR  AU.  fPTTFRTA  GO  FTRf  ANOTHFR  ROMND, 

C 

TFfWiwi  .GT.O,  )GnTnR7 
TFfPASSA'.GT.O.  IGOTODT 
C 

C  TF  THIS  RDljOn  JS  A  KILL  FOR  ML  fWTTFPTA  GO  TO  MF^T  SA^^PL^. 

C 

GOTORO 

are  used  to  transfer  control  to  fire  another  round  or  to  proceed  to  run 
the  next  sample  case  after  examining  the  kill  criteria  and  target  passen¬ 
ger  flags. 

The  statements 


c 

c  TF  THIS  TS  THF  'aAIN  POUND  GO  TO 
C 

1  TFT  THio'.NF.OTGOTOPOS 
C 

C  TF  THF  SPOTTTNG  PTFLF  HAS  SHOT  MOPF  OP  FQUAL  TO  THF  5PFCTFTF0  POUNDS. 
C 

TFT  THI!  .GF.NROSI  TGOTOPOT* 

C 

c  WAS  THTS  POUND  SFNSFD? 

C  . 

TFfPANDOMfn.GT.PP  SIIGOTOPno 
C 

C  rOPRFfT  FT  TN  PAPALl FLOGP AM  FASHTON, 

C 

TALI.  NRAN^t  (PI  *  PP) 
xr 1 aXr 1 -X 1 ♦XC-Pl *STGXS1 
VC  1 =YC 1 -V 1 f Yr-RP*ST GYS 1 
FACsO. 

C 

C  FTRF  ANOTHFR  SPOTTING  POUND. 

C 

GOTOPOn 


are  used  to  identify  the  round  fired.  If  the  round  was  a  main  round, 
control  is  passed  to  Statement  203.  If  the  round  was  a  spotting  rifle 
round,  a  test  for  the  maximum  number  of  spotting  rounds  fired  is  made; 
then  a  uniform  random  number  is  chosen  to  determine  if  the  miss  of  the 
spotting  rifle  round  was  sensed.  If  the  miss  of  the  spotting  rifle 
round  was  sensed.  Subroutine  NRAN31  is  called  to  furnish  the  random 
normal  deviates  required  to  update  the  horizontal  and  vertical  coordi¬ 
nates  of  the  center  of  impact  of  the  spotting  rounds.  Control  is  then 
transferred  to  Statement  200  to  fire  another  spotting  round. 
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MAIN  ROUTINE  (10  of  13) 
The  statements 


c 

C  CnRRFCT  THF  rFMTPR  OF  TARGFT  OUT  OF  fFNTFR  OF  IMPACT  BFCAUSF  IT  WILL 

C‘  BF  AODFH  TM  WHFN  MAIM  POHMO  IS  FT»FD. 

C 

?ou  xnsxn-vr 

YCtsYrt-vr 

noTopoi 

are  used  to  modify  the  horizontal  and  vertical  coordinates  of  the  center 
of  impact  for  the  spotting  rounds  before  control  is  transferred  to 
Statement  201  to  fire  a  main  round. 

The  statements 


?0h  TMl?rTHi;>4^t 

THI ^STHI ^♦THI 1 

c 

C  TF  ALL  OF  SPnTTTWG  RTFl.F  RO  IMHS  HAVF  BFF*^  FTRFH  AND  ROUNn  IS  NOT  WITH 
C  IN  1?  TNCHFS  OF  TARGET  GO  TO  NFXT  SAMplF  OTHERWISE  FIRE  MAIN  ROUND. 

C 

TFfAXfl  l*t?..GT.yn  GOTORR 
TF(AyfNXl  +  l?..lT.Xn  GOTOQR 
TF{AVfNYl  +  t?..LT.Yn  GHTORO 
TFfAYfl  )-lR..GT.Yn  GOTORQ 
GOTORfia 


are  used  to  update  the  spotting  round  counters  and  to  determine  whether 
to  proceed  to  run  the  next  sample  case  or  to  fire  a  main  round  if  the 
spotting  round  impact  point  lies  within  12  inches  of  the  enveloping 
target  grid  cell  area. 

The  statements 


c 

C  IF  IMF  MISS  nF  THE  MAIM  RntJMD  IS  NOT  SENSED  FTRF  ANOTHER  ROUND  WITH 
C  NO  rORRFfTIOM, 

C 

IFfRANDOMfl 1 .GT.PRORS)  GO  TO  Q? 

C 

C  CORRECT  CT  FOR  NEXT  ROUND  USING  PARALLELOGRAM  MpyHOO. 

C 

CALL  npam^I (RI #9?) 
xri=xri-xi4xr-Rt*STGxs 
YC1=Yri-Yt^VC-RP*STr,YS 
FACsO, 

R7  rONTTMUF 


are  used  to  update  the  horizontal  and  vertical  coordinates  of  the  center 
of  impact  for  the  main  rounds  after  calling  Subroutine  NRAN31  to  provide 
the  required  random  normal  computational  deviates  (in  accordance  with 
the  equations  which  follow)  if  the  miss  of  the  main  round  was  sensed. 
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MAIN  ROUTINE  (11  of  13) 


(2-5) 


(2-6) 


The  statements 


rr  f  .0*  iGOTOo^ 

c 

C  rHFTK  FII7F  ACTMATIOM. 

C 

TFf  TH?0,FO.O)GnTn^O;> 

c 

C  PFWSONMFL  rASMAUTTFS  CHMWlJTFn. 

C 

SMMCAS  sSlIMCA-S  ♦rPK(K??.S)/PASSM)  aCPASSM-SUMCAS  ) 
^0?  7r j.s)37f j.si^suvrAS 
Qft  rONTTNUF 

c 

c  FTWF  NFXT  main  POUND. 

c 

TF( J.L T.NRDSlGOTOPflS 
C 

C  GO  TO  NFXT  SAMPLF. 

C 

9Q  rONTTNlIF 


are  used  to  check  the  passenger  flag  in  order  to  compute  the  number  of 
target  passengers  killed  for  the  cell  hit  under  consideration  if  the 
fuze  actuation  flag  is  not  equal  to  zero.  The  cumulative  Pk  array  con¬ 
taining  target  passengers  killed  per  rounds  fired  is  then  updated  before 
either  firing  another  main  round  or  continuing  with  the  execution  of 
another  Monte  Carlo  trial. 

The  statements 


c 

c 

c 

7 

C 

c 

c 


sampling  TS  COMPLFTF. 

'^PTTFrfe.  n 

FnWMATr//lPH  NPDS  PK/) 

COMPMTF  MFAN  PK  AND  PH  TAPGFT  DATA 


AND  /JRTTF  RFSULTS. 
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MAIN  ROUTINE  (12  of  13) 

SAWP-msamP 
SAMPTsI ,/samp 

on  Q6  T?51 .mrhs 

0n91 T=1 »h 

91  7rT-?.T)  =  7f  T?,T)*SAMPT 

«7  WPT  TF(^,  #»)!>•(  7(  !?•  T  ).  Txl 
9h  rnwTTNHF 

are  used  to  reduce  the  cumulative  Pk  target  array  representing  9,604 
Monte  Carlo  samples  for  all  kill  criteria  considered  to  mean  values  and 
to  write  them  for  output.  The  average  number  of  spotting  rounds  fired 
and  the  number  of  times  that  all  the  spotting  rounds  were  fired  is  also 
written  by  the  execution  of  the  statements 


IHt  ^StHl  ^/NS»MP 

c 

C  TM1?  TS  NUMHFP  OF  TTMFS  F-IPFO  ALL  SPOTTIMG  RnUN^S*  IHM  IS  AVFRAGF 
C  MMMSFP  OF  SPOTTTMG  R'lU'ViHS  FTPFH. 

C 

WRTTF(6,?Sl  )  TH1?. TH1 ^ 

FOPMATrpTSI 


The  statements 


c 

C  rnv4PnTF  MFAN  PK  ANH  PH  TT'^F  OATA  AMD  WPTTF  RESULTS* 

c 

WRT  TF(h.#i1  ) 

nnisairi  pf>i 
nnisajri,f> 

ts«  AKTl  (  T  ,,nr  AKTLf  I  #  !)*SAMPl 

nniss  1x1 

wRTTF  f8. 1  f  AKTl  f  I  ».n  .  Txl  ,ht  ) 

FORMAT ( 1 OFfl.S) 

WPTTFf#,.1S7) 

1S7  FOR^^TflMO) 

1SS  rnMTTMilF 


are  used  to  reduce  the  cumulative  Pk  time  array  to  mean  values  and  to 
write  them  for  output  in  accordance  with  equation: 
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^  r>  o  r>  ri 


MAIN  ROUTINE  (13  of  13) 
The  statements 


tO'^PMTF  AVFPAGF  OF  Pns.  FTPFP. 

TFf  MOPTI  .FO*^)r,OTr|P 
CAU  KTIPAT 

r,n  TO  RFAH  Npyr  data  sft. 
nnrn? 

h  FOPMATfTS. 

81  FOPMATflMI)  • 

H?  FORMATf/Z^IH  M»IMRFR  OF  SA^PLFS  2  ,T6./lHn 
FND 

are  used  to  call  Subroutine  KILRAT  to  compute  the  average  number  of  main 
rounds  fired  per  case  if  the  program  option  flag  is  not  equal  to  three. 
Control  of  the  program  transfers  to  Statement  2  to  read  a  new  data  set. 
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FUNCTION  CNF(X) 


This  function  is  used  to  compute  the  cumulative  normal  distribution 
values  needed  to  determine  the  probability  of  target  cell  impact  by  a 
main  round  for  the  closed  form  solution. 

The  statements 


FIIMCTION  CNF(K) 
piiRpnsr 

rn^^OLATTV^  normal  DTSTPTRUTTIN  VAlUFS  FOR  CFLL  PH  HSFD 
tM  rinSFO  FDRM  SOLUTION. 


c 

F  =  0, 

AXsAPSf  X)  *  ,7071 
TFf  Ax.r,F.s. 
F=.S-.S«FRF(AX) 

I  TFf  X.r,F.0.)Fs1  .-F 

TNFsF 
RFTURN 
FND 


are  used  to  generate  values  from  the  cumulative  normal  distribution 
utilizing  the  error  function  associated  with  the  normal  curve.  Control 
then  is  returned  to  the  calling  routine. 


oon  nno  non  nnr>  nor» 


SUBROUTINE  INPUT  (1  of  5) 


The  purpose  of  this  subroutine  is  to  read  and  write  the  input  data, 
to  convert  the  input  data  into  units  suitable  for  program  execution,  and 
to  set  up  the  target  grid  and  Pk  arrays.  The  subroutine  is  begun  by  the 
statements 


•SUB^ntiTTNF  T'VtPllT 


C  PURPOSF. 

C  RFAO  AMO  i^RTTE  INPUT  DATA,  CONVFRT  UNITS  OF  INPUT  DATA,  AND  SFT  UP 
C  THF  TAROFT  fiPTF)  AND  PK  ARRAVS. 

C 

nTMFNSTON  PKt^OnO  ,KI,AX(Pnn,AY(Sl  I  ,C  ASF(aOI  ,  Z  f  tP#P)  »Z1  f  61  ) 
COMMON  PK. AX, AY,SX»,SYR,SXB,SYW,SXS,SYS,PRS#XR,YR,XC, Yr#RFLf TFl , 

1 TFS, XM1 ,XM?. YnFF.NOPTl , MROS , NX . N Y , FL T , S T t ,ST?, AMTl , AMT ? , T K (SO # ?00 ) 
?. 7. Z1 , AKTL(61 ,91 .STRXL.STnYL.RFLF 

■^..S  XSI.S  YSl.S  XRI.S  YRl.S  XRl.S  YR 1  ,  PRS  !  ,  XR  I  ,  YR 1  ,  FL  T  t  ,  TF  It  , 
^TFS1,XMtt,XMPl,AMT1i,AMTP1,STtl,ST^1, SIGXL t ,STGYl 1 .NRDSl .PASSN 


which  are  used  to  allocate  storage  for  arrays  and  to  provide  exchange  of 
data  between  this  subroutine  and  the  remainder  of  the  program. 

The  statements 


RFAO  AND  WRTTF  TTTLF  CARD. 

RFAOfS.PO.FNDsfe)  CASF 

?0  FORMATfunAP} 

WRTTF(6,1R)  CASF 
1«  format  flHI  .iJOApi 

RFAD  main  ROUND  DATA. 

RFAOrs, 1 ,FNn=6)NOPTl , NRDS • MOP T . XT , YC , SXS , SYS , SXR , SYR , SXR , SYR , PRS , R 
IFl  , YDFF. XR. YR.FLT, TF1 , TF S , XM i , XM? , AMT  I , AMT? , ST  1 , ST?, STGXL , STGYL # 
?RFl F, PASSN 

CHFCK  FLAG  TO  RFAD  SPOTTING  ROUND  DATA. 

TF  (NOPT  1  ‘.NF.?IGOTn?0  I 

RFADfS, POO ) NRDSl .SXS1,SYSt.SXR1,SYRl,SXRl,SYRl,PRSt,XRl,YH!,FLT1, 

1 TF1 1 . TFSt , XMt 1 , XM?1 , AMTIt , AMTPt ,STn • ST?t , STGXL 1 , STGYL 1 
?00  format r TS.SX, 7F 1 O.P, /RFI 0. P, /RF 1 0 .P) 

NUMRFR  or  MATN  ROUNDS  FIRFO  PFR  SAMPlF  SPFCTFTFD, 

?01  TF(NRnS^I.  T,  1  0)  NRDS=10 
TFfNRDS.GT. 1 9)  NR0S=19 
TFfNOPTI  .FQ. NRnSn 

WRTTF  MATN  ROUND  INPUT  DATA. 

WRTTF(6,?1N0PT1 ,NRDS,NnPT 

?  FORMATf  ///////////PI  H  ST^^UI.ATTON  s  ,  T  S  , /?0  M  Ml  |MRf  r  OF  W 

10UNOS  =  ,TS,/«H  NOPT  =  ,TSt 
WRTTF(6, %) 
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SUBROUTINE  INPUT  (2  of  5) 

I  FOR^ATf HORT/nviTAl.  V^RTTCAL  ,//.^1X,7m  MFTFRf^,  Sy,7H  MF 
t  TFR5i,/) 

wPTTFfh,U)SXS,SYS.SXH.SY«.SXR,SYR,  V  n  ,  YH  .  X  C  .  Yf  ,  5  T  G  XL  .  S  T  T,  YL 
a  FORWATf^KH  SFNST'jr,  FRRORS  (  S  T  r>  nrv)  ,  \  ?  .a , /?iiH  VARTARl.F  RIAS  (ST 

10  OFv)  •  ?x»?F  1  ?.a. /?aH  RAMrvn^  frwhrs  (Sto  orvi  ftxfo 

?  HTAS,  1SX,?F1P.4|./11  M  ATM  CFNTFs  .  1  SX,  ?F  1  /,  1  OH  lAY  FRWnR,lf»X, 

^PFIP.U/) 

WRTTF  ((»,7)PRS.RFL,  YDFF.RFLF 

7  FORMATfiRH  PRORARTLTTY  HF  SFNSTNO  MISSING  RIUMO  s  ,  F  1  ? ,  6  , /P’^H  RFLIA 
1RTLTTY  OF  ROUNO  s  .FIP.^./^iJH  VFRTTCAL  OFFTLAHFn  r^^R^M^'ATFS  s  • 
?F1P.^».7H  TNCHFS/.P5H  RFL  T  ART  1. 1  T  Y  OF  FIJ7F  r  .FPB.6) 

WRTTF(B,1RIFLT,TF1 ,TFS#XM1 ,xmp,aMT1 ,AMT?,sT1 ,ST? 

19  FnPMATr/////^?X. 1BH  TTMF  0ATA(SFCS),///18H  TTMF  OF  FLIGHT  =  ,F1?.6 
1.///P8X,inH  1ST  R0tlN0.RX.16H  SORSFOUFMT  R0S..//2SH  FIXED  MEDIAN  TT 
?MF  .F1P.6.10X.F1P.6./PSH  VARIARLE  MEDIAN  TJMF  ,F1?.6.10X. 

■^F1?,6./PSH  MINIMUM  TIME  .  F  1  ?  .  6  ,  t  0  X  ,  F  1  ?  .  6 , /?SH  SLOPE  OF 

aOTSTRIROTTON  . F 1 ? . 6 . 1 0 X . F 1 P . 6 1 

are  used  to  set  the  number  of  main  rounds  to  be  fired  per  sample  and  to 
read  and  write  the  input  data  for  the  main  rounds  and,  if  the  program 
option  flag  so  governs,  for  the  spotting  rifle  round  data. 

The  statements 


c 

C  convert  MATM  PD.  DFI.TVFRY  data  from  METRIC  TO  ENGLISH  STANDARD  UNITS. 
C 

SXSsSXSaXR. ^7 
SYS=:SYS*X0.'^7 
SxR=SXRaXR.X7 
SVRsSVR*  XR. ^7 
SXRzSXR*XO. X7 
SYRzSYR*  XR. ^7 
STGXLrSTGXl  *XR.'i7 
STGYl sSTGYL*  ^R.^7 
XR  zXR  aXR,'^; 

YR  zYR  XQ.  M 

xcsxr *XR. 57 
YCsYr*5R.57 
C 

C  CHECK  FLAG  FOR  i^RTTTNG  SPOTTTNG  RTFLF  ROUND  DATA. 

C 

IFINOPTl .NE.PIGOTOPnp 
RFLFlzl. 

RFL 1  =  1  . 

wRTTF{6.1Ri 

ARTTF(6.Pn51 

PnX  FORMAT! ///?0H  SPOTTING  RIFLE  DATA, //I 
.NRTTF(6.5) 

WRT  TF(6.RI SXSl .SYSl .SXW1 .SYRl , SXRt .SYRl . XRI , YRl . X C . VC , S I GXL 1  * S T GYL 
11 

WRTTF(6. 7IPPS1 .RFL1 . YOFF.RFIFl 

•YPTTF  (6,  IRIFLTl  .  TFl  I  ,  TFSl  .  XMt  1  ,  XMPI  .  AMJ  U  •  AM  T  p  t  ,  S  T  11  .STPl 

C 

C  COMVFPT  SPOT  PD.  DFLTVFRY  DATA  FROM  METRIC  TO  ENGLISH  STANDARD  UNITS. 
C 

SXSlzSXSl*^R.57 
8YS1=SYS1*^R, ^7 
SXR1zSXRM5R.57 
SYRlrSYRl *59.57 
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SUBROUTINE  INPUT  (3  of  5) 

SXR1=SXP1 *^9.^7 
SYR1 rSYWI * ^7 
srr.xi  1  =sTr*xLi  *  jq.  ^7 
STHYL  1  =STr,YH 
xRi  =xm  ♦ 

YMIsYRl S7 

are  used  to  convert  the  input  data  from  metric  units  to  English  standard 
measure  units.  The  spotting  rifle  round  data  is  written  if  the  option 
flag  is  equal  to  two. 

The  statements 


c 

C  rWFW  PpnFTrTFNfY  Ty^F  FACTORS  SFT. 

c 

TFfSTl  1  .FQ.ft.l  Sn 
TFfSTPI  .FQ.O.)  ST21s.4iqrt^ 

20P  TF(ST1  .FQ.O.  ) 

TFf  STP.FO.n.  )  STPr.^lQS^ 

C 

C  CHFrK  TNPilT  C^)^^TROL  FLAG  FOP  IJSTMR  PPFVTOMS  HATA  SFT  HP. 

C  .  . 

TF  fNnPT.FQ.O^PFTUPN 


are  used  to  set  the  slope  of  the  variable  time  distribution  for  the  crew 
proficiency  exponential  equation  and  to  return  control  to  the  calling 
routine  if  the  flag  NOPT  is  equal  to  zero  indicating  that  the  vulnera¬ 
bility  matrix,  the  identifying  kill  array,  and  the  Pk  array  remain  as 
in  the  previous  input  data  set. 

Descriptive  target  data  such  as  the  number  and  size  of  target  grid 
cells  in  the  horizontal  and  vertical  directions  and  the  minimum  horizontal 
and  minimum  vertical  coordinate  values  which  can  be  assumed  by  the  grid 
area  enveloping  the  target  are  read  by  the  execution  of  the  statements 

c 

C  RFAH  TARr.FT  DATA. 

C 

RFAPf  S,  1  >NX,^^Y 

MXl=NX+i 

MYIsNV+t 

1  FORMAT  f  ^TS,SX,6F1 O.P, /8F I  0 . P , /«F 1 0 . ? , /8F 10.?) 

IS  RFAn(S.SnXMT»y.  YMTN,DX.nv 


The  statements 
c 

C  SFT  VUL NFRARTLTTY  CARO  ITMtT. 
C 

NCs'^noo 

c 

C  SFT  IIP  VIM.MFRARTUI  TY  MATRIX, 

C 

AXn  )rXMTM 
on  18  TsP.MXI 
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SUBROUTINE  INPUT  (4  of  5) 

\h  AX(T)rA^f  T-l 

ft  V  r  n  sYMTM 
nn  17  T=P.NV1 
17  AVnirftYf  T-1  )+DY 

c 

C  tNTTTALT7F  TD  KILL  ARRAY. 

C 

nnS«Tst,NX 
nOSajsl.MY 
sa  TKf.J.T)  =  0 

are  used  to  limit  the  number  of  Pk  input  cards  to  3,000,  to  initialize 
the  identifying  kill  array,  and  to  set  up  the  arrays  containing  the 
incremented  horizontal  and  vertical  values  of  the  grid  area  enveloping 
the  target. 

The  statements 


c 

C  FNTFR  VUl NFRARTLTTY  CARR  LOOP. 

C 

nnsPTsi.wr 

Mri=s 

TFfPASSN‘.I.F.0.1Mri=a 

c 

C  rMFTK  FLAG  FOR  STANOARH  ARRAY  SFT  UP. 

c 

TFfMnPT.FO.l IRDTO^O^ 

C 

C  RFVFRSF  TARGFT  nWTFNTATTHM. 

C 

RFAOfS.  ^Oa>  X,  Y,  (PK  (  T  .  f  1  ..1=1  •MCt  1 

Yr-X 

GHTD^O? 

50*^  RF  AO  rS.  ^OSIX.Y.  (PK  (  T  .J)  .1=1  .Nr  !  "I 
mu  FORMAT  (PF7.P.  ^RX  ,  MF  7  .  n 
50S  FnRMAT(PF7.P.  8X.SF7,n 

c 

C  CHFTK  FOR  LAST  VHNF  RARTl  T  T  Y  CAR^. 

C 

m?  TF  (X.GT.RRR.  )  RFTIJRN 

C 

C  C TNVFRT  TARGFT  COOROTNATFS  INTO  TO  KTLI  ARRAY  SURStRTPTS. 

C 

.f1s((X-XMTN>/nX)  +  1  . 

J?=((Y«YMTN)/OY)f1 . 

SP  TKflP.jnsT 
in  FORMAiriHi)  • 

RFTIIRN 

Si  FOPMATraFIO.! .TS) 

#»  rALl  FXTT 
FNO 


are  used  to  read  the  Pk  vulnerability  data  for  the  various  selected  kill 
criteria  for  a  given  point  on  the  target.  If  the  personnel  casualty  flag 
is  less  than  or  equal  to  zero,  only  four  target  kill  criteria  are  con¬ 
sidered.  Otherwise,  personnel  casualty  data  are  also  read.  If  the 
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SUBROUTINE  INPUT  (5  of  5) 

input  control  flag  is  equal  to  one,  the  Pk  and  identifying  kill  arrays 
are  set  up  in  standard  fashion  so  that  the  given  target  coordinates  are 
converted  into  the  identifying  kill  array  subscripts  which  correspond 
with  the  subscripts  of  the  target-enveloping  grid  area  arrays.  If  the 
input  control  flag  is  greater  than  one,  the  orientation  of  the  target  is 
reversed,  and  the  appropriate  Pk  vulnerability  data  are  read  accordingly. 
Control  is  then  returned  to  the  calling  routine. 
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SUBROUTINE  KILRAT  (  1  of  2) 


This  subroutine  is  used  to  compute  and  write  the  average  number 
of  main  rounds  fired  for  each  case. 


The  statements 


SUHRfUJTTNF  KTL»Ar 


C 


«ft****A*******«A«*** 


C  PMffPnSF 

C  ClWPNTfS  AWn  WRTTFS  THF  AVFRAGF  NUMRFR  OF  MAIN  ROUMDS  FTRFD  PFR  CASF 

. . 


OTMFNSTON  AKTLf  8I#P>»F(  6n#S^^K(IP  *P)fAA(  6 1  )  #  PK  f  0  0  #  ^  1  « 

iAxfPot).AYrsn 

rnMMfiN  PK.AX,  AY,STGXR.STr,YP,STnXH,  S!GYB,SinxS,STGYS,Sir,  •XB^YB, 

1  xr  ,  YC  .PFl  •  TF1  ,  TFS.  XVI  ,  X*^? .  YOFP .  NOP  T!  ,NR  ,  NX  ,  MY  ,  FL  T  ,  vST  t  ,  ST?,  AMT  1  , 
?AMT?,TKfSO,?On)  ,SM»<,F,  AKTl  ,  S  T  GX  L  ,  S 1  G  Y  L  .  RFL  F 

^i.STGXSi  .STGYS1  .  S  T G X B  I ,  S  T  G Y R  1  ,  S  t  GXR1  ,  STGYH  1  ,  PRSl  ,  XB  i  ,  YB  1  ,  FI  T  !  .  TF  1  1  . 
iiTFSl  ,XMi  1  .  ,  AviTl  1  .  AMT?1  .ST  n  .ST?1  ,STGXL  )  ,STGYL  1  .NROSl  .PASSN 


are  used  to  allocate  storage  for  arrays  and  to  provide  for  the  exchange 
of  data  between  this  subroutine  and  the  remainder  of  the  program. 

The  statements 


FNTFR  r ASF  I nnp. 

nn  ?  T 1 =1 .6 
SMTKrSMK (1 .  I  n 
AMlrP. 

AMRsNR 

IrNR 


FNTFR  RDIINO  LOOP  AND  COMpglE  AVG.  RDS.  FTRFO, 

nnpsTsP.J 

SMTKsSMTKf ( AUL* ( SMK ( T . It ) -SMK ( T*1  ,  T  1  n  ) 

?S  AtiLsAilL*!. 

OTG  ALL  SAMPlFS  RFGTSTFR  A  KTLL? 

TFfSMKfNR.  T1  KGT..PR9Pt  GOTOGP 


are  used  to  enter  the  case  loop  and  to  execute  the  round  loop  to  compute 
the  average  number  of  main  rounds  fired  for  a  given  case  or  kill  cri¬ 
terion  in  accordance  with  the  equation: 


r 

AVI  =  2  i(Pi  -  Pa  ■,) 
i=l  ^ 


(2-13) 
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SUBROUTINE  KILRAT  (2  of  2) 

If  a  kill  occurs  within  the  set  number  of  rounds  fired  for  each  sample 
of  a  case,  control  transfers  to  Statement  99  to  write  out  the  results. 

The  statements 

A»sS,- (smk  ( j-s,  I  n  f  ,  1 1 1  ( j- t  1 ( j-?»  t  1 1  ^smk  ( j-i  #111) 

11  * 

XXsSmkCJ  .Tn-SMK(J-S#Tt  1 
TF(xx,uF.n.)r,fUO<)«) 

A/sAR/XX 

c  COMPUTF  AVG  Rns  FIPFn  TMClunT^JG  NOM-K  T  LL  T  MG  SAMPLFS  AMD  FX  TRAPOL  AT  TOM  FACTOR 
C 

AVr(A74-ANR)*(1  .1  •  TH  1  ♦SMTK 

GOTnQR 
C 

c  WRTTF  AVG  ROUMOS  FTRFO  PFR  CASr. 

C 

RR  AVsS'^TK 
RR  'jgPTTFCR.PRl  AV.  T1 

FOR^^ATflRH  AVG.  MH.  ROS  s  •F7.?,8H  CASF  =  •7S.//1 
?  rOMTTMilF 
WF  TURN 
FNO 


are  used  to  compute  the  average  number  of  rounds  fired  per  case  including 
those  rounds  fired  on  non-killing  samples  and  an  error  estimate  based 
upon  extrapolation  to  the  limit  in  accordance  with  the  equation: 


(2-14) 


After  all  results  are  computed  and  written,  control  is  returned  to  the 
calling  routine. 
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SUBROUTINE  NRAN31  (XI,  X2) 


The  purpose  of  this  subroutine  is  to  generate  random  numbers 
representing  normal  deviates  for  use  in  establishing  aimpoints,  impact 
points,  and  time-to-kill  values  for  both  main  and  spotting  rounds. 

The  statements 


SIJUBO'JTTNF  VRAM^KXt.*?) 

PURPIISF 

RFMFRATFS  PANDH'^  NORMAL  nFVIATFS. 

FIINCTTON  RANDOM  RFQUT^FD. 

**«**A**ft4A*A««**«««***********«4*A*A********«*****************«<l**'***********^ 

)(^  =  SORT(-?.^ALOG(*<ANOOW(t  ))) 

=  8S‘%0  7?*«AMnnM(  1) 

XP  =  X  ^♦STNf 
XI sx  ^*rns( xai 
RFTURM 
f  Njr> 


are  used  to  produce  two  normal  deviates  which  are  employed  to  establish 
aimpoints,  impact  points,  and  time-to-kill  values  after  calling  Function 
RANDOM  and  applying  the  sine  and  cosine  functions  as  factors.  Control 
is  then  returned  to  the  calling  routine. 
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FUNCTION  RANDOM  (N) 


This  function  is  used  to  generate  uniform  random  numbers  in  the 
interval  between  zero  and  one.  The  subroutine  is  begun  by  the  statements 

FIIMCTTOM 

C*********A*******************A<»***»****<^************<k*****^***************^*^** 

C  PuwpnsF 

C  GF^^FRATFS  U^^TFnRM  PAVnOV  NUN.iRFRS  RFTwFFM  0.  AMO  1. 

C*<k^>*****»****<k***A-A**<k************A****AA  ************************************** 

c 

rriMMON/  PAM^M/tX 


which  provide  for  the  exchange  of  data  between  this  function  and  the 
remainder  of  the  program. 

The  statements 


TXr  IX*hSS  %<i 
rfdx'is.h.h 

s  T^  =  TX  +  Pl/i7a8^6  47  +  t 

6  sTX 

1  SUMrStJM*‘.u8S#»81  ^F-R 
PA^fOnM  =  S^M 
RFTnWM 
FMD 


generate  a  random  number  by  first  multiplying  two  large  integers.  If 
overflow  is  produced,  the  random  number  is  made  positive  by  adding 

Ol 

2  .  The  resulting  random  number  produced  after  subsequently  multiplying 

oi 

by  the  constant  2  ,  lies  in  the  interval  between  zero  and  one.  Control 

is  then  returned  to  the  calling  routine. 
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SUBROUTINE  SSHOT  (PCH,  PCHl )  (1  of  3) 


The  purpose  of  this  subroutine  is  to  compute  the  probability  of 
hitting  the  target  for  a  single  main  round  and  the  probability  of  hitting 
the  target  for  a  single  main  round  given  a  sensed  miss  by  using  a  closed 
form  solution. 

The  subroutine  is  begun  by  the  statements 

- . SUBPOUTTMF  sswm (PCH, PCHl ) 

c 

c  PHPPnsP 

C  PM  TARGFT  A  STNGLF  ROU’^O  AND  FOP  A  ST^»GLF  MAIN  ROUND 

C  GTVFN  A  SFNSFD  MISS  RY  USING  A  Cl.nSFO  FORM  SOLUTION, 

C 

C  FUNGTinN  CNF  RFQUIPFD, 

C 

;  OTMFNSTON  PKf^OOO  ,  8  I  •  A  Y  I  ;>0  n  ,  A  Y  (  S  1  I  *  7  (  1  P »  9  )  #  7  t  f  6 1  ) 

rOMMON  PK. AX, AY^STCyPtSIGYR, STGXR,  S T GYR , S T GX S , S T GYS , PROPS • XB , YR, 

1 xr . YCtRFL. TF t • TFS, XMt , XM? , yOFF , NOPT I , NROS, NX • N Y , FL T • ST  1 , ST?, AMT  1 , 
PAMTP.TKfSO  .?00),7,71.AKTI  f 61 , 9), STGXL . 5 TGYL , RFLF 

STGXS1 ,STGYS1 , STGXRl , S I G YR 1  . S T G X R 1 • S I G Y R 1 , PR S t , X R 1 , YR 1 , FL  T  t , T F  1  1  , 

^TFSI . XMn ,XM?1 , AMT ll,AMT?l.ST11.5T?t,S!GXLt, STGYL 1 ,NROSl ,PASSN 


which  are  used  to  allocate  storage  for  arrays  and  to  provide  for  the 
exchange  of  data  between  this  subroutine  and  the  remainder  of  the 
program. 

The  standard  deviations  of  impact  error  in  the  horizontal  and 
vertical  directions  are  computed  in  accordance  with  the  equations: 


II 

X 

D 

/  2  .2  .2 
rXR 

II 

>- 

/  2  .2  .2 
fayR  +  c^yB  '^YL 

i - 1 

(2-17) 


(2-18) 


(2-19) 
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SUBROUTINE  SSHOT  (2  of  3) 


(2-20) 


These  standard  deviations  of  impact  error  utilize  the  system  and  random 
input  delivery  error  data  for  a  main  round  and  for  a  main  round  given  a 
sensed  miss  by  execution  of  the  following  statements 


PCM  sO. 

PCHIsO, 

COMPIITF  IMP&CT  SICfMAS  FT?  A  ST^r.I.F  '^ATM  STtMO. 

str.K=soRT(STr.xp**?  +  sir.xR**?  +  siRXL**?) 
Str.YsSQRT  fSirtVP**?*STGYR**?  +  S1GYL**?) 
SIRAsS0RT(SIGXS*»?y?.*.S10X«**?) 

STGHrSOPT (SIGYS**?  +  ?.*StGYH*  *?) 


The  statements 

CHFCK  FOP  SPOTTIMG  ROUND  DATA. 

IFfNnPTI .NF.2)GnTO^ 

MODIFY  impact  sigmas. 

SIGXsSDRT  fSIGX**?Y( .S*(SIGXS1 ySIGXL  t ) ) *  * P  +  SIGXRl ** ?) 
SIGV=SQRT(SIGY**?K .S*(SIGYS1 ♦RTGYL1 n**?+SIGVRl**?) 


are  used  to  modify  the  standard  deviations  of  the  main  round  computed 
above  by  incorporating  spotting  rifle  round  delivery  error  data  if  the 
option  flag,  NOPTl ,  is  equal  to  two  in  accordance  with  the  equations: 


Q 

X 

li 

2  2  2 
ay  +  [•5(a)^5^  +  Oj^li)]  +  <^XR1 

Q 

-< 

II 

'2  2  2 

Qy  ■*"  ^  ^YRl 

(2-21) 


(2-22) 
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SUBROUTINE  SSHOT  (3  of  3) 

The  statements 

FNTFR  l.nnp  Tn  SFLFCT  CFLL  VFRTTCAL  rOOPDlNATFS, 


I  OnpJsKNY 

YPsAYf J>1 )-YC 
Y1=AYCJ  )-Yr 

TMFCK  TARRFT  nFFlLAOFD  f  ni^O  T  T  T  O^^S  . 

TFf  AYf.14.1  ).LF.YnFF)  GOTO? 

!  TFfAYfJ  ).IT.Y0FF)  Vl=vnFF-Yr 

FNTFR  l.nnP  TO  SFLFCT  CFLL  HnRT/HMTAL  mORO  T  WA  TFS , 

001  1  =  1 

CHFCK  THAT  CFLL  TARGFT. 

TFf  TKf.|.T)*LF.O  ir,OTni 
X1=AX(T  )-XC 
YP  =  AXfT  +  n-XC 

COMPUTF  ANO  SUM  CFLL?  =>H*S. 

PA  =(rNF((XP-YR)/STGV)-CNF ( ( X 1  - YB 1 / S T G X ) 1  * ( C NF ( ( Y ?- Y R WS T G Y ) -C NF ( ( 

1 YI-YB)/STGY1 ) 

PA1 =f CNF  f X?/STGA1-CMF (XI /STGA ) ) * (CMF ( Y^/STGRI -C NF ( Yl /SIGRl ) 

PCM  rPCH  fPA 
PfHI sPCHI ♦PAI 
t  rOMTTMirF 
?  C^^JTTNUF 
RFTURN 
FNO 

are  used  to  execute  the  target  cell  loops  in  order  to  set  up  the  vertical 
and  horizontal  coordinates  of  each  target  cell  to  be  examined.  If  the 
target  cell  is  defiladed,  the  vertical  cell  loop  is  re-entered.  If  the 
target  cell  has  no  associated  Pk  values,  the  horizontal  cell  loop  is 
re-entered.  The  probabilities  of  hitting  each  cell  for  a  single  main 
round  and  the  probabilities  of  hitting  each  cell  for  a  single  main  round 
given  a  sensed  miss  are  computed  by  using  the  cumulative  normal  distri¬ 
bution  function,  then  appropriately  summed  in  accordance  with  Equations 
2-15  and  2-16  before  control  returns  to  the  calling  routine. 


(2-15) 


a 

Ph  =  S 
i=l 

PA. 

(2-16) 
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ABBREVIATIONS  AND  SYMBOLS 

The  following  pages  are  a  dictionary  providing  definitions  for  the 
various  abbreviations  and  symbols  employed  in  the  Simulation  Model.  The 
list  is  ordered  as  follows.  Those  symbols  appearing  in  blank  COMMON 
appear  first,  in  alphabetical  order.  These  are  followed  by  those  symbols 
appearing  in  labeled  COMMON.  The  labeled  COMMON  groups  are  ordered  alpha¬ 
betically  by  COMMON  label  name  and  those  symbols  within  each  COMMON  block 
are  ordered  alphabetically-  Finally,  those  symbols  which  appear  only  in 
a  single  routine  are  grouped  according  to  the  routine  in  which  they  appear. 
The  MAIN  Routine  appears  first  with  the  subroutine  groups  ordered  alpha¬ 
betically  according  to  subroutine  name.  The  symbols  within  each  subrou¬ 
tine  group  are  ordered  alphabetically. 

Along  with  the  symbol  definition,  the  equivalent  symbol  in  the 
Mathematical  Model  is  provided  for  cross  referencing.  The  parameter 
units  are  provided  as  well*  A  dash  (---)  in  the  units  column  indicates 
that  the  parameter  is  nondimensional . 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (1  of  7) 


Equivalent  in 


Abbreviation 
or  Symbol 

Mathematical 

Model 

Definition 

Units 

AKIL(61 ,9) 

Array  containing  mean  time 

Pk  values  and  mean  time  Ph 
values  (cumulative  values 
until  reduced  to  mean  values). 

AMTl 

— 

Minimum  time  for  first  main 
round. 

seconds 

AMT2 

— 

Minimum  time  for  main  rounds 
subsequent  to  the  first  round. 

seconds 

AMTll 

— 

Minimum  time  for  first 
spotting  round. 

seconds 

AMT21 

— 

Minimum  time  for  spotting 
rounds  subsequent  to  the 
first  round. 

seconds 

AX (201) 

— 

Array  containing  the  incre¬ 
mented  X  values  of  the  grid 
area  enveloping  the  target. 

inches 

AY(51) 

“  “  ” 

Array  containing  the  incre¬ 
mented  Y  values  of  the  grid 
area  enveloping  the  target. 

inches 

FLT 

Projectile  time  of  flight  of 
main  round. 

seconds 

FLTl 

— 

Projectile  time  of  flight  of 
spotting  round. 

seconds 

I K( 50, 200) 

Identifying  kill  array  with 
subscripts  matching  AX  and  AY 
arrays  containing  sequence 
numbers  associating  Pk 
values  with  target  cells. 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (2  of  7) 

Equivalent  in 
Abbreviation  Mathematical 

or  Symbol  Model 

NOPTl 

NR 

NRDS 

NRDSl 

NX 

NY 

PASSN 

PK{3000,8) 

PROBS  — 

PRS 

PRSl 


Definition 

Flag  governing  program 
options: 

=1 :  main  rounds  only; 

=2:  spotting  rifle  data 
included; 

=3:  number  of  rounds  fired 
per  sample  is  one. 

Number  of  main  rounds  to  be 
fired  per  sample. 

Number  of  main  rounds  to  be 
fired  per  sample. 

Number  of  spotting  rifle 
rounds  to  be  fired  per  sample. 

Number  of  target  cells  in 
horizontal  (X)  direction. 

Number  of  target  cells  in 
vertical  (Y)  direction. 

Number  of  target  passengers. 

Array  containing  probability 
of  target  kill  input  values 
for  designated  x-  and  y- 
coordi nates. 

Probability  of  sensing  miss 
of  main  round. 

Probability  of  sensing  miss 
of  main  round. 

Probability  of  sensing  miss 
of  spotting  round. 
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Units 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (3  of  7) 


Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

REL 

— 

RELF 

— 

SIGXB 

°XB 

SIGXBl 

— 

SIGXL 

°XL 

SIGXLl 

°XL1 

SIGXR 

X 

D 

SIGXRl 

°XR1 

SIGXS 

°XS 

SIGXSl 

°XS1 

SIGYB 

"^YB 

Definition 

Reliability  of  main  round. 

Reliability  of  fuze  of  main 
round. 

Standard  deviation  of  hori¬ 
zontal  variable  bias  for  main 
rounds. 

Standard  deviation  of  hori¬ 
zontal  variable  bias  for 
spotting  rounds. 

Standard  deviation  of  hori¬ 
zontal  lay  error  for  main 
rounds. 

Standard  deviation  of  hori¬ 
zontal  lay  error  for 
spotting  rounds. 

Standard  deviation  of  hori¬ 
zontal  random  error  for  main 
rounds. 

Standard  deviation  of  hori¬ 
zontal  random  error  for 
spotting  rounds. 

Standard  deviation  of  hori¬ 
zontal  sensing  error  for  main 
rounds. 

Standard  deviation  of  hori¬ 
zontal  sensing  error  for 
spotting  rounds. 

Standard  deviation  of  vertical 
variable  bias  for  main  rounds. 


Units 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 


inches 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (4  of  7) 


Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

Definition 

Units 

SIGYBl 

— 

Standard  deviation  of  vertical 
bias  for  spotting  rounds. 

inches 

SI6YL 

°YL 

Standard  deviation  of  vertical 
lay  error  for  main  rounds. 

inches 

SIGYLl 

°YL1 

Standard  deviation  of  vertical 
lay  error  for  spotting 
rounds. 

inches 

SIGYR 

"^YR 

Standard  deviation  of  vertical 
random  error  for  main  rounds. 

inches 

SIGYRl 

°YR1 

Standard  deviation  of  vertical 
random  error  for  spotting 
rounds. 

inches 

SIGYS 

"^YS 

Standard  deviation  of  vertical 
sensing  error  for  main  rounds. 

inches 

SIGYS1 

"^YSl 

Standard  deviation  of  vertical 
sensing  error  for  spotting 
rounds. 

inches 

SMK(19,9) 

Array  containing  mean  Pk 
values  for  selected  kill 
criteria  and  mean  Ph  of 
target  for  rounds  fired. 

STl 

ST 

Slope  of  time  distribution 
for  first  main  round. 

— 

ST2 

ST 

Slope  of  time  distribution 
for  main  rounds  subsequent 
to  the  first  round. 

STll 

ST 

Slope  of  time  distribution 
for  first  spotting  round. 

— 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (5  of  7) 


Equivalent  in 


Abbreviation 
or  Symbol 

Mathematical 

Model 

Definition 

Units 

ST21 

ST 

Slope  of  time  distribution 
for  spotting  rounds  subse¬ 
quent  to  the  first  round. 

“  “  — 

SXB 

^XB 

Standard  deviation  of  hori¬ 
zontal  variable  bias  for 
main  rounds. 

inches 

SXBl 

--- 

Standard  deviation  of  hori¬ 
zontal  variable  bias  for 
spotting  rounds. 

inches 

SXR 

^XR 

Standard  deviation  of  hori¬ 
zontal  random  error  for  main 
rounds. 

inches 

SXRl 

— 

Standard  deviation  of  hori¬ 
zontal  random  error  for 
spotting  rounds. 

inches 

SXS 

""xs 

Standard  deviation  of  hori¬ 
zontal  sensing  error  for 
main  rounds. 

inches 

SXSl 

— 

Standard  deviation  of  hori¬ 
zontal  sensing  error  for 
spotting  rounds. 

inches 

SYB 

°YB 

Standard  deviation  of  vertical 
variable  bias  for  main  rounds. 

inches 

SYBl 

— 

Standard  deviation  of  vertical 
variable  bias  for  spotting 
rounds. 

inches 

SYR 

'^YR 

Standard  deviation  of  vertical 
random  error  for  main  rounds. 

inches 

SYRl 

— 

Standard  deviation  of  vertical 
random  error  for  spotting 
rounds. 

inches 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (6  of  7) 


Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

Definition 

SYS 

^YS 

Standard  deviation  of  vertical 
sensing  error  for  main  rounds. 

SYS1 

Standard  deviation  of  vertical 
sensing  error  for  spotting 
rounds. 

TFS 

TFS 

Fixed  median  time  for  main 
rounds  subsequent  to  first 
round. 

TFS1 

TFS 

Fixed  median  time  for  spotting 
rounds  subsequent  to  first 
round. 

TF1 

— 

Fixed  median  time  for  first 
main  round. 

TFll 

— 

Fixed  median  time  for  first 
spotting  round. 

XB 

XB 

Horizontal  fixed  bias  for 
main  rounds. 

XBl 

— 

Horizontal  fixed  bias  for 
spotting  rounds. 

XC 

XC 

Intended  horizontal  aim 
center  of  target. 

XMl 

XM 

Variable  median  time  for 
first  main  round. 

XM2 

XM 

Variable  median  time  for  main 
rounds  subsequent  to  the 
first  round. 

XMll 

XM 

Variable  median  time  for 

f i rst  spotti ng  round . 
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Units 

inches 

inches 

seconds 

seconds 

seconds 

seconds 

inches 

inches 

inches 

seconds 

seconds 

seconds 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


COMMON  (7  of  7) 

Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

Definition 

Units 

XM21 

XM 

Variable  median  time  for 
spotting  rounds  subsequent  to 
the  first  round. 

seconds 

YB 

YB 

Vertical  fixed  bias  for  main 
rounds. 

inches 

YBl 

— 

Vertical  fixed  bias  for 
spotting  rounds. 

inches 

YC 

YC 

Intended  vertical  aim  center 
of  target. 

inches 

YDEF 

— 

Vertical  defiladed  coordi¬ 
nate. 

inches 

Z(19.9) 

— 

Array  containing  mean  Pk 
values  for  selected  kill  cri¬ 
teria  and  mean  Ph  of  target 
values  for  rounds  fired. 

mm  mm  ^ 
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COMMON/RANDM 

Abbreviation 
or  Symbol 

II,  IX 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


Equivalent  in 
Mathematical 

Model  Definition 

Random  seed,  initially  equal 
to  11111111. 


Units 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


MAIN  ROUTINE  (1  of  3) 


Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

C(4) 

--- 

FAC 

— 

HWM 

HWM 

HWl 

— 

HW2 

— 

IHW 

— 

IHIO 

— 

IHll 

— 

IH12 

— 

IH13 

— 

IH20 

— 

J 

— 

JIO 

Definition 

Array  containing  flags  indi¬ 
cating  kills  for  up  to  four 
established  kill  criteria. 

Flag  controlling  incorpora¬ 
tion  of  lay  error. 

Incremental  time  factor  for 
variable  crew  proficiency. 

Flag  indicating  achievement 
of  all  kill  criteria. 

Uniform  random  number  used 
in  determining  target  kills. 

Flag  indicating  first  main 
round  hit  of  each  sample. 

Flag  indicating  main  round 
being  fired. 

Spotting  rifle  round  counter 
for  each  sample. 

Number  of  times  all  spotting 
rounds  fired. 

Average  number  of  spotting 
rounds  fired. 

Flag  indicating  fuze  actua¬ 
tion  on  round  hit. 

Counter  for  main  rounds  fired 
per  sample. 

Intermediate  value  identify¬ 
ing  main  round  fired. 


Units 


seconds 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


MAIN  ROUTINE  (2  of  3) 

Equivalent  in 
Abbreviation  Mathematical 

or  Symbol 

Model 

Definition 

Units 

JIO 

--- 

Intermediate  first  time  Pk 
element  value  in  AKIL  array. 

seconds 

K22 

Relevant  Pk  values  for 
the  corresponding  target 
cell. 

NSAMP 

— 

Program  sample  size. 

— 

PCH 

Probability  of  hit  of  single 
round  (closed  solution, 

Monte  Carlo  not  used). 

... 

PCHl 

Probability  of  hit  for  single 
round  given  a  sensed  miss 
(closed  solution). 

^ 

R1 

—  —  — 

Random  normal  deviate  used  in 
computing  impact  points  of 
main  and  spotting  rounds. 

R2 

— 

Random  normal  deviate  used  in 
computing  impact  points  of 
main  and  spotting  rounds. 

^  ^ 

R3 

■■■ 

Random  normal  deviate  used  in 
computing  impact  points  of 
spotting  rounds. 

"  “ 

R4 

Random  normal  deviate  used  in 
computing  impact  points  of 
spotting  rounds. 

SAMP 

— 

Program  sample  size. 

— 

SAMP  I 

— 

Inverse  of  SAMP. 

— 

SUMCAS 

— 

Number  of  target  passengers 
killed. 

— 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


MAIN  ROUTINE  (3  of  3) 


Abbreviation 
or  Symbol 

Equivalent  in 
Mathematical 
Model 

Definition 

Units 

T 

Time  of  round  impact. 

seconds 

T1 

Random  normal  deviate  used  in 
updating  time  for  main  and 
spotting  rounds. 

—  —  “ 

T2 

Random  normal  deviate  used  in 
updating  time  for  main  and 
spotting  rounds. 

XCl 

XCl 

Horizontal  center  of  impact. 

inches 

XI 

XI 

Horizontal  impact  point  of 
main  cr  spotting  rounds. 

inches 

X2 

Horizontal  coordinate  of  lay 
error  contribution. 

inches 

YC1 

YCl 

Vertical  center  of  impact  for 
main  or  spotting  rounds. 

inches 

Y1 

Y1 

Vertical  impact  point  of 
main  or  spotting  rounds. 

inches 

Y2 

— 

Vertical  coordinate  of  lay 
error  contribution. 

inches 

4-39 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


FUNCTION  CNF 


Abbreviation 

Equivalent  in 
Mathematical 

or  Symbol 

Model 

Definition 

AX 

— 

Absolute  value  of  X. 

F 

— 

Represents  value  of  cumulative 
normal  distribution. 

X 

— 

Value  for  which  cumulative 
normal  distribution  is  solved. 

LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  INPUT  (1  of  2) 

Equivalent  in 
Abbreviation  Mathematical 
or  Symbol  Model 

CASE 

DX 

DY 

01 

02 

NC  — 

NCI 

NOPT 


Definition 
Title  of  run. 

Target  cell  size  in  the  hori¬ 
zontal  direction  (X). 

Target  cell  size  in  the  ver¬ 
tical  direction  (Y). 

Index  of  IK  array  matching  X 
value  for  input  Pk's  with  the 
index  of  the  AX  array  in  the 
horizontal  direction  for  a 
given  target  cell. 

Index  of  IK  array  matching  Y 
value  for  input  Pk's  with  the 
index  of  the  AY  array  in  the 
vertical  direction  for  a  given 
target  cell . 

Maximum  number  of  input  Pk 
cards  which  can  be  read 
(3,000). 

Represents  the  number  of  Pk 
values  per  card  to  oe  read  for 
a  given  x,y  target  coordi¬ 
nate. 

Input  control  flag. 

NOPT  =  0:  AX,  AY,  IK,  and 
Pk  arrays  remain  as  in  pre¬ 
vious  input  data  set; 

NOPT  =  1 :  above  arrays  set 
up  in  standard  fashion; 

NOPT  ^2:  reverses  standard 
target  orientation. 

Number  of  target  cells  plus 
one  in  the  horizontal  (X) 
direction. 


Units 

i nches 

inches 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  INPUT  (2  of  2) 

Equivalent  in 
Abbreviation  Mathematical 

or  Symbol  Model 

NYl 

RELl 

RELFl 

X 

XMIN 

Y 

YMIN 


Definition 

Number  of  target  cells  plus 
one  in  the  vertical  (Y) 
direction. 

Spotting  round  reliability. 

Spotting  round  fuze  reliabi¬ 
lity. 

Represents  horizontal  (X) 
target  coordinate. 

Represents  minimum  horizontal 
value  of  grid  area  enveloping 
target. 

Represents  vertical  (Y)  target 
coordinate. 

Represents  minimum  vertical 
value  of  grid  area  enveloping 
target. 


Units 


inches 

inches 

inches 

inches 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  KILRAT 


Equivalent  in 


Abbreviation 
or  Symbol 

Mathematical 

Model 

Definition 

ANR 

r 

Number  of  main  rounds  fired 
per  sample. 

AR 

... 

Intermediate  variable  used  in 
computing  average  number  of 
rounds  fired. 

AUL 

i 

Intermediate  factor  repre¬ 
senting  the  number  of  rounds 
fired  in  computing  SMIK. 

AV 

AVI 

AV2 

Average  number  of  main  rounds 
fired  to  achieve  a  kill  for 
a  designated  case. 

AZ 

M  M  M 

Intermediate  variable  used  in 
computing  average  number  of 
rounds  fired. 

11 

— 

Case  number. 

J 

r 

Number  of  main  rounds  to  be 
fired  per  sample. 

SMIK 

AVI 

Average  number  of  main  rounds 
fired  to  achieve  a  kill  for  a 
designated  case  excluding 
non-killing  trials. 

XX 

— 

Intermediate  variable  used  in 
computing  average  number  of 
rounds  fired. 

LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  NRAN31 


Equivalent  in 


Abbreviation 
or  Symbol 

Mathematical 

Model 

Definition 

XI 

— 

Random  normal  deviate. 

X2 

— 

Random  normal  deviate. 

X3 

— 

Normal  random  number. 

X4 

— 

Uniform  random  number, 
0  <  X4  <  2Tr 

Units 


radians 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


FUNCTION  RANDOM 


Equivalent  in 

Abbreviation 

Mathematical 

or  Symbol 

Model 

Definition 

Units 

SUM 

— 

Uniform  random  number, 

0  <  SUM  <  1 . 

— 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  SSHOT  (1  of  2) 


Equivalent  in 


Abbreviation 
or  Symbol 

Mathematical 

Model 

Definition 

Units 

PA 

PAi 

Probability  of  hitting  a 
given  grid  cell  of  target 
with  a  single  round. 

PAl 

PA. 

Probability  of  hitting  a 
given  grid  cell  of  target 
with  a  single  round  given  a 
sensed  miss. 

PCH 

Ph 

Probability  of  hitting  target 
with  a  single  round. 

— 

PCHl 

Ph 

Probability  of  hitting  target 
with  a  single  round  given  a 
sensed  miss. 

SIGA 

Horizontal  standard  deviation 
of  impact  error  for  a  single 
main  round  given  a  sensed 
miss. 

inches 

SIGB 

°B 

Vertical  standard  deviation 
of  impact  error  for  a  single 
main  round  given  a  sensed 
miss. 

inches 

SIGX 

^X 

Horizontal  standard  deviation 
of  impact  error  for  a  single 
main  round. 

inches 

SIGY 

^Y  ' 

Vertical  standard  deviation 
of  impact  error  for  a  single 
main  round. 

inches 

XI 

— 

Represents  lower  horizontal 
(X)  value  of  target  grid  cell 
being  examined. 

inches 

X2 

— 

Represents  greater  horizontal 
(X)  value  of  target  grid  cell 
being  examined. 

inches 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
(SIMULATION  MODEL) 


SUBROUTINE  SSHOT 

Abbreviation 
or  Symbol 

Y1 

Y2 


(2  of  2) 

Equivalent  in 
Mathematical 

Model  Definition  Units 

—  Represents  lower  vertical  (Y)  inches 

value  of  target  grid  cell 
being  examined. 

—  Represents  upper  vertical  (Y)  inches 

value  of  target  grid  cell 
being  examined. 


Next  page  is  blank. 
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SECTION  V 


INPUT 


The  input  data  card  set  contains  information  concerning  target 
description,  target  vulnerability,  delivery  accuracy,  round  character¬ 
istics,  and  crew  proficiency. 

The  paragraphs  that  follow  provide  general  information  regarding 
input  data  preparation.  An  illustration  of  the  input  data  deck  organi¬ 
zation  (Figure  5-1)  and  descriptions  of  the  input  data  card  formats  are 
provided  at  the  end  of  this  section. 

The  program  option  flag,  NOPTl ,  controls  execution  of  three 
possible  alternatives: 

•  When  NOPTl  =  1,  only  main  rounds  are  fired  for  each  trial  of 
the  simulation  program. 

•  When  NOPTl  =  2,  main  and  spotting  rifle  rounds  are  fired  for 
each  trial  of  the  simulation  program. 

•  When  NOPTl  =  3,  only  one  main  round  is  fired  for  each  trial 
of  the  simulation  program. 

The  input  control  flag,  NOPT,  also  governs  three  options: 

•  When  NOPT  =  0,  the  vulnerability  matrix,  the  identifying  kill 
array,  and  the  Pk  array  remain  as  constituted  in  the  previous 
input  data  set. 

•  When  NOPT  =  1,  the  vulnerability  matrix,  the  identifying  kill 
array,  and  the  Pk  array  are  set  up  in  the  standard  fashion 

to  read  the  first  group  of  kill  data  from  a  vulnerability  card. 

•  When  NOPT  ^  2,  the  sign  of  horizontal  target  coordinates  is 
reversed,  thus  reversing  standard  target  orientation.  This 
option  causes  a  second  group  of  kill  data  to  be  read  from  a 
vulnerability  card. 

The  NRDS  flag  controls  the  maximum  number  of  main  rounds  to  be 
fired  for  each  sample: 

•  If  the  input  value  of  NRDS  <  10,  a  maximum  of  10  main  rounds 
will  be  fired  per  sample  except  that  the  program  option  flag, 
NOPTl,  prevails  for  one  main  round. 

•  If  the  input  value  of  NRDS  >  19,  a  maximum  of  19  main  rounds 
will  be  fired  per  sample. 
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Warhead  input  also  includes  round  and  fuze  reliability  data. 

The  target  is  defined  by  delineating  its  maximum  and  minimum  grid 
cells  in  both  the  horizontal  and  vertical  directions.  The  number  of 
horizontal  and  vertical  target  grid  cells  and  the  size  of  these  cells 
complete  the  input  data  description  of  the  target.  The  number  of  vul¬ 
nerability  input  cards  is  limited  to  3,000.  The  last  of  these  is 
denoted  by  a  value  exceeding  999.  For  a  given  point  on  the  target, 

Pk  cell  data  for  each  designated  target  kill  criterion  are  read  and 
stored  for  later  processing.  The  program  accommodates  up  to  four  target 
kill  criteria  and  a  fifth  category  concerning  target  passenger  kills  for 
an  APC  type  target. 

Main  round  and,  if  applicable,  spotting  round  delivery  input 
parameters  for  a  stipulated  range  include  intended  aimpoint,  standard 
deviation  of  horizontal  and  vertical  sensing  errors,  standard  deviation 
of  horizontal  and  vertical  variable  bias,  standard  deviation  of  hori¬ 
zontal  and  vertical  random  errors,  standard  deviation  of  horizontal  and 
vertical  lay  errors,  horizontal  and  vertical  fixed  bias,  and  the  vertical 
defilading  coordinate  if  the  target  is  defiladed.  All  of  the  above 
delivery  input  parameters  are  quantified  in  metric  units  except  the  ver¬ 
tical  defilading  coordinate  which  is  described  in  English  standard  units. 
Time  factor  parameters  for  both  main  and  spotting  rifle  rounds  include 
projectile  time  of  flight,  fixed  median  times,  variable  median  times, 
minimum  times,  and  slope  of  time  distribution  data. 
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SECTION  VI 


SAMPLE  CASE  OUTPUT 


This  section  contains  a  sample  problem  with  a  data  checksheet 
(Figure  6-1)  to  illustrate  the  input  variables  which  are  printed  out 
in  Figure  6-2  in  the  same  order  as  they  appear  in  Figure  6-1.  Although 
there  are  over  1,000  vulnerability  cards  required  for  this  particular 
sample  case,  only  a  representative  few  are  shown  in  Figure  6-2.  (The 
entire  vulnerability  card  deck  appears  in  Appendix  A.)  This  section 
also  provides  samples  of  the  output  data  formats. 
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DATA  CHECKSHEET 


Figure  6-1.  Sample  Case  Input  Data  Sheet 
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Figure  6-2.  Sample  Case  Card  Input  Print-out 
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PRINTED  OUTPUT 


The  printed  data  include  various  input  data  parameters;  closed 
form  solutions  for  Ph  single  round  data;  mean  target  Pk,  mean  target  Ph, 
mean  time  Pk,  and  mean  time  Ph  data;  and  the  average  number  of  rounds 
fired  per  case. 

Input  Data  Listing 

The  data  listed  in  Figure  6-3  represents  input  data  parameters 
including  the  run  title,  program  and  input  option  flag  indicators, 
number  of  main  rounds  fired  per  sample,  the  system  and  random  delivery 
errors  in  the  horizontal  and  vertical  directions,  the  probability  of  a 
missing  round  being  sensed,  round  reliability,  the  vertical  defilading 
coordinate,  and  fuze  reliability.  Time  data  including  projectile  time 
of  flight  and  crew  proficiency  time  factors  are  followed  by  the  number 
of  Monte  Carlo  trials  run  and  the  results  of  the  closed  form  solution  to 
the  probabilities  of  hitting  the  target  for  a  single  main  round  and  for 
a  single  main  round  given  a  sensed  miss. 

Target  Results 

The  first  10  rows  of  column  1  in  Figure  6-4  indicate  the  number  of 
rounds  fired  per  Monte  Carlo  trial.  The  next  four  columns  of  10  rows 
represent  the  mean  target  Pk  values  for  the  designated  kill  criteria. 

Mean  target  Ph  data  are  presented  in  column  7  in  which  the  first  answer 
for  one  round  fired  per  sample  should  show  agreement  with  the  results 
of  the  closed  form  solution.  Row  11  of  columns  1  and  2  depict  the  number 
of  times  throughout  the  samples  that  all  the  specified  maximum  number 
of  spotting  rifle  rounds  were  fired  and  the  average  number  of  spotting 
rounds  fired  per  trial. 

For  those  APC  target  data  sets  including  passenger  kill  information, 
column  5  is  used  to  print  the  results  of  the  average  number  of  personnel 
casualties  per  rounds  fired. 
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range  ■  250  MEIERS  NON  DEFILADED 


SIMULATION  B  1 

number  of  rounds  ■  10 

NOPT  ■  1 


HORIZONTAL  vertical 
METERS  meters 


0,0 

0.13'#T 
0.320n 
-0.4430 
0.0 
0.0 

probability  of  sensing  missing  round  =  0.0 

reliability  of  hound  =  1.000000 

vertical  defiladed  coordinates  =  -999.000000  inches 

reliability  of  TiJZE  =  1.000000 


sensing  errors  (STD  DEV)  0,0 

variable  bias  (STD  DEV)  0.1240 

Random  errors  (Std  DEV)  6.32on 

FIXED  bias  -0,5293 

AtM  center 

LAV  ERROR  0*0 


TIME  OaTA(SECS) 


TIME  OF  flight 


2.000000 


1ST  ROUND 


subsequent  ROS 


FIXED  MEDIAN  time  0.0  0.0 

Variable  median  time  15.000000  15.000000 

MINIMUM  time  0.0  0.0 

SLOPE  OF  distribution  0.0  ('.0 

ppoe.  OF  HIT  OF  FIRST  ROUND  =  0.9954607 

PROS.  OF  HiT/SENSED  MISS  =  0.9464230 


number  of  samples  ■  9hn4 


Figure  6-3.  Sample  Case  Input  Data  Print-out 
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TARGET  KILL  CRITERIA 


NO.  1 

NO.  2 

NO.  3 

NO.  4 

Target 

NROS 

•  Ph 

Pk3 

. ^^4  — 

— 

-  --  Ph 

1 

0.7229279 

0.6665972 

0.7338608 

0*5867347 

0.0 

0.9958349 

Z 

0.9203456 

0.««32777 

0.9246147 

0.0278841 

0.0 

0,9999999 

3 

0.9801124 

•  0.9593918 

0.9818825 

0*9282590 

0.0 

0.9999969 

4 

0.0930236 

0.9819866 

0.9936484 

0,9668887 

0.0 

0,9999999 

5 

0.0977092 

0.9917741 

0.9978133 

0,9844856 

0.0 

0.9999999 

6 

0.0991f,69 

0.9960432 

-  ■0.999<i710 

"0.9919824 

T) ,  0  " 

0.9999999 

7 

0.09966^6 

0.9984381 

0.9996876 

0.9969803 

0.0 

0.9999999 

6 

0.9997917 

0.9993752 

0.9997917 

0,9990628 

0.0 

0.9999999 

9 

0.9996958 

0.9906876 

0.99989^38 

0,9994793 

0.0 

0.9999999 

10 

0.0998958 

0.9996876 

0.99989^58 

0,9995834 

0.0 

0.9999999 

j-►o  o-«-(Average  number  of  .  spotting  rifTe.  rounds  fired_p.er_tri^ 

(Number  of  times  throughout  the  samples  that  all  the  specified  number 
of  spotting  rifle  rounds  were  fired) 


Figure  6-4.  Mean  Pk  Results 
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Time  Results 


For  the  given  kill  criteria,  each  of  the  elements  of  the  time  Pk 
arrays  listed  in  Figure  6-5  represents  a  2-second  interval  from  time 
0  to  120  seconds  reading  from  left  to  right.  The  mean  time  Pk  arrays 
for  the  various  kill  criteria  are  ordered  sequentially  with  the  first 
array  denoting  the  first  kill  criterion,  etc.  However,  the  sixth  array 
is  associated  with  column  7  of  the  target  results  and  pertains  to  mean 
Ph  time  data. 

Round  Expenditure 

The  final  output  results  concerning  the  average  number  of  main 
rounds  fired  to  achieve  a  kill  per  case  are  shown  in  Figure  6-6.  The 
last  case  represents  the  average  number  of  rounds  to  merely  hit  the 
target. 
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Time  Pk  Array  for  Kill  Criterion  1  (0-120  Seconds) 

0.0  0.0  0.0  0.00e6U  0.02405  0. 07486  0. 14031  0.23719  0.32309  0.40254 

0.47053  0.53405  0.59298  0.64202  0.68347  0.71855  0.75531  0.78488  0.80789  0,83090 

0.84044  0.86724  0.88286  0.89733  0.90816  0.91681  0.92774  0.93419  0.94242  0.95033. 

0.95637  0.96147  0.96533  0.96855  0.97220  0.97470  0.97709  0.97959  0.98188  0,98407 

0.98573  0.9869b  0.9887b  0.99052  0.99167  0.99209  0.99344  0.99406  0,99448  0.99531 

0.99583  0.99625  0.99677  0.99688  0.99750  0.99781  0.99813  0.99833  0.99B54  0,99865 

0,99885 

Time  Pk  Array  for  Kill  Criterion  2  (0-120  Seconds) 

0.0  0.0  0.0  0.00187  0,02166  0.06799  0.13723  0.21949  0,29988  0.37328 

0.43669  0.49573  0.55144  0.59B50  0.63921  0.67368  0.71127  0.74136  0.76510  0.78946 

0.81018  0.82861  0.84527  0.86131  0.87453  0.88505  0.89617  0.90650  0.91660  0.92420 

0,93180  0.93857  0.94419  0.94877  0.95325  0,95731  0.96210  0.96564  0.96887  0.97189 

0.97459  0.97678  0.97918  0.98126  0.98386  0.90428  0.98563  0.98667  0.90803  0.98928 

0.99063  0.99105  0.99209  0.99292  0.99334  0.99386  0.99427  0.99459  0.99500  0.99594 

0.99625 

Time  Pk  Array  for  Kill  Criterion  3  (0-120  Seconds) 

0,0  0.0  0.0  0.00260  0.02436  0.07580  0.15160  0.24063  0.32799  0.40848 

0.47720  0.54102  0.59985  0.64890  0.69023  0.72511  0.76145  0,79134  0.81331  0.83611 

0.85454  0.87182  0.88713  0.90108  0.91139  0.91982  0.93076  0.93711  0.94544  0.95262 

0.95866  0.96366  0.9b74l  0.97032  0.97376  0.97626  0.97865  0.98095  0.98313  0.98532 

0,98698  0.98823  0.98990  0.99177  0.99292  0.99323  0.99406  0.99469  0,99511  0.99594 

0.99646  0.99667  0.99719  0.99719  0.99761  0.99792  0.99823  0.99844  0.99865  0.99875 

0.99896 

Time  Pk  Array  for  Kill  Criterion  4  (0-120  Seconds) 

0.0  0,0  0.0  0.00156  0.01853  0.06029  0.i2099"0.19336  0.26395  0,32997 

0.38776  0.44138  0,49386  0.5387J  0.87820  0.61329  0,65025  0,68201  0./0929  0,73449 

0.75625  0.77707  0.79394  0,81195  0.82799  0.84059  0.85579  0.86662  0.87953  0.88848 

0.89838  0.90764  0.91524  0.92149  0,92763  0.93326  0.93888  0.94461  0.95002  0,95398 

0.95793  0.96137  0.96481  0.96751  0.97137  0.97251  0.97543  0.97678  0,97855  0.98022 

0.98209  0.9835b  0.98501  0.98636  0.98740  0.98803  0.9890^0.9896^  0.99084  0.99209 
0,99261  " 

0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

0,0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

0,0  0.0  0.0  0,0  0.0  0.0  0.0  _ 0,0^  0,0  0,0 

0.0  0.0  -  0.0  0.0 .  0.0  0.0  0.0  0.0  '  0.0  0.0 

0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0,0  0,0. 

0,0  0.0  0.0  0,0  0.0  0.0  0.0  0,0  0,6  .0,0 
0.0 

Time  Ph  Array  (0-120  Seconds) 


0.0 -  0.0  . '  0.0  "  0.00375  0.03374  0.10235  0.20429  0.32507-0.44086  0.54894 

0.63515  0.71064  0.77468  0.82257  0.86141  0.09035  0.91451  0.93308  6.94648  0.95845 
0.96699  0.97386  0.98063  0.98396  0.98698  0.98896  0.99115  0,99188  6.99344  0.99459 
0.99563  0.99625  0.99688  0,99740  0.99771  0.99792  0.99844  0,99375  6,99885  C. 99917 
0.99927  0.99938  0,99969  0,99969  0.99979  0.99979  0,99979  0.99979  0,99979  0.99979 
0.99979  0.09979  0.99990  0.99990  0.99990  0.99990  0.99990  0,99990  0.99990  0.99990 
0.99990 


Figure  6-5.  Mean  Time  Pk  Results 
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Kill  Criteria 


AVG.  NO. 

Rr>s  « 

1*39 

Case  = 

-  1 

av6*  no* 

ROS^  « 

i.w 

Case  = 

2 

'  - tl0~  2  - 

AVG.  NO. 

ROS  * 

1.37 

Case  b 

3 

NO.  3 

AVG.  NO. 

ROS  « 

1.72 

Case  = 

4  . . 

-  mr  4-  - 

AVG.  NO. 

RDS  a 

0.0 

Case  ■ 

5 

— 

AVG*  NO* 

RDS  m 

1.00 

Case  « 

‘6  - 

- ?h- 

Figure 

6-6. 

Case  Rounds  Fired 

Next  page  is  blank. 
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SECTION  VII 


SOURCE  LISTING 


The  following  section 
Routine,  Subroutine  INPUT, 
Subroutine  SSHOT  presented 


contains  the  source  listings  for  the  MAIN 
Subroutine  KILRAT,  Function  RANDOM,  and 
in  Figures  7-1  through  7-7. 
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onn  non  nnnn  nnn  non  nnnnnnnnnnnnnnnnnnnn 


.IwrM/SS  MDnrt  FOR  niRFCT  FTRF 

MAIN  PPnnRAM 
PMRPnSF 

inCATFS  TMPACT  P0TNT3  f>F  MAIN  AMO  SPOTT!^^G  RTFLF  ROUNDSJ  MnWITDRS 

THF  TIMF  TO  KTI.LI  COMPilTFS  PK  AMO  PH  TARGET  DATA  AND  PK  AND  PH  TIME  DATA 


SlIRoPnGPAMS  RFOIJTPFD 
INPUT-  RFAOS  INPUT  DATA. 

KIIRAT-  rOMPUTFS  AND  w/RTTFS  THF  AVERAGE  NU'^RFR  OF  MAIN  RDUMDS  FIRED 
PFP  TARF. 

NRAM^t-  PROVIDFS  NHRn^AL  DFVTATFR  FOR  MONTE  CARLO  SAMPLING 
RANDOM-UNTFfIRM  RANDOM  NUMRFR  GENERATOR. 

SSHOI-  rOMPIlIFS  PH  FO^  A  SINGLE  maTN  ROUND  AND  PH  FDR  A  SINGLE  ROUND 
GIVEN  A  SENSED  MISS  USING  A  CLOSED  FORM  SOLUTION. 


nTMFNSTON  PKf'^OOO  »  AX  f  ?0n  ,  AYfSl  IfZCIP  ,  P  )  #  C  ( ft  )  •  Z  1  C  1  ) 

rOMMON  PK. AX. Ay,STGXR,STGYR,STGXfl,  S T G Y B  ,  S T G X S . S I G YS .  PRORS  .  X B  .  Vfl . 
ixr.  YC.RFL.  TF1  .TFS.XM1  ,  XM? ,  ¥  DF  F  .  •'•OPT  I .  MRDS .  N  X  .  N  Y  .  FL  T  .  S  T  1  .  S  T2  .  AM  T  1  . 
;>AMTp,TK(Sn.?nn),7,71  .AKILC^I.R)  .STGXL  ,  STGYL  ,RELF 

^.STGXSI .S7GYS1 .STGXBl .STGYR1  .SIGXRl . SIGYRl ,PRS1 .XRl , YBt .FLTt , TFl 1 . 
iJTFS1.XMH.XM?i.AWTt1.AMT?1.STn.ST?1,  ST  GXLl  .  STGYL  I  .NRDSl  .  PASSN 

rOMMON/  RANDM/Tl 

SFT  Hr  TO  ANY  ODD  NUMRFR  FOR  RANDOM  SFFD, 

I i St  1  1  1  11  1  ! 

RFAD  INPUT 
?  CALL  INPUT 

COMPOTE  PROBABILITY  OF  HIT  OF  FIRST  ROUND  AND  PROBABILITY  OF  HIT  GIVEN 
A  SENSED  MISS  IN  A  CLOSED  FORM  SOLUTION  AND  WRITE  RESULTS. 

r A1  I  SSHOT (PCH.PCHl ) 

WRTTFfh.ftO)  PCH.PCHl 

fcO  FORMATf/HH  PROR.  OF  HIT  nf  FIRST  ROUND  r  ,  F  1  8 . 7 .  /  ?« PROR  .  OF  HIT 
l/SFNSFD  MISS  =  ,Flft.7) 

initialize  sample  ST7F. 

NSAMps9f,oa 

INITIAL T7r  TIMF  ARRAY  A NO  PPOR A R T »  I T Y  ARRAY. 

DOPSJrl .ft 

noisftin.fti 

ISft  AKit  fT..nrO. 

DO  PS  Trt .NRDS 


Figure  7-1.  Source  Listing,  MAIN  Routine  (1  of  7) 


c;iJuo  ouu  fcjuu  .ouuuuu  ouo  cjuouuo 


<>S  7rT.J)20. 

PRTMT  nilT  SAVPLF  3T7F. 

W9TTF(#».#»?)M3AMP 

start  SAMPI  TNG 

TMt?sO 
TH! S=0 

m  99  Tsi .NSAMP 

TNTTTALT7F  VARTARLFS  FOR  FACH  SAMPLF. 

SIIMCAS  =n. 
yrizo, 
vriso. 

THlOsO 
TH1 IrO 
T=n. 

THj^sO 
nn9a.isi.R 
9a 

TF  nMLV  maTN  rounds  CONSIDFRFO  GO  TO  ?01. 

TFT  Nnprr.NF.PTGOTIPOT 
CALL  NRANSl (Rt .R?) 

'SFT  T  TO  TTMF  Of  FLIGHT  OF  ONF  SPOTTING  RTFLF  SHOT. 

TsFI  T1 

•  XCI.YCI  TS  CFNTFR  OF  IMPACT  FOR  FIRST  ROUND  OF  SPOTTING  RIF|  F  SHOTS. 

yC1 sWI ^STGXBl ♦XCfXRI 
Yri=R?*STGYBlfYC+VR1 
FAr=1 . 

DRA./^  NORMA!  DFVIATFS  FOR  RANDOM  AND  LAY  FRRORS  OF  SPOTTING  RTFLF. 

20n  CALL  NRAN^I fRl ,R?) 

CALL  NRAN^I (R^,Ra) 

THIt  TDFNTTFTFS  THF  SPOTTING  RTFIF  ROUNDS, 

TH11=TH1i+1 

DRAW  NnRMAI  DEVTATFS  "OR  TT‘^F  FOR  SPHTTING  RIFLF. 


r 

c 

c 

c 


CALL  NRAN-^t  (T1  ,  T?) 

FDR  FTPS!  ROUND  TI^F  OF  SPOTTING  PTFIF  GO  If)  ?10. 
TFf  TH1 1  .FD. 1 IGOTOPI 0 


Figure  7-1.  Source  Listing,  MAIN  Routine  (2  of  7) 
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r»or>  ooo  oooo  noo  ooo  r»r»r»  orvo  nno  rir>r>  r»or> 


C  POR  SUftSFQMr^T  ROUND  TTMFS, 

C 

HWMsywpl  *F^(P(T  1  *.ST^n 
TFrnWM j  r,AV!T?l  )  HW'^sAMT?! 

TsTfHAM^TFSl  , 

ooTnpi 1 

?10  1  *rxp  (T  1  AST  t  n 

TF  f  HinM.I  F  ,  A'^TI  n  HAlMsAMTIt 

TrTfH»^M  +  TF1  1 

rnMTT'^UF 

TyPATT  PniNT  FOR  THIS  SPOTTTNG  RTFI.F  SHOT. 

XI sXC  t  aRI aSTGXRI aR^^STOXL  t  aFAC 
Y1  =  Yn  ♦RPaSTGVWI  fRUASIGYM  aFAC 

GO  SFF  TF  THTS  SHOT  HIT. 

GnTni#,o 

ORA-V  RANDOM  OFVTATFS  FOR  VARIARLF  BIAS  ERRORS  FOR  FTRST  MAIN  POUND. 
?01  TALL  NRAN^I f Rl .RP) 

UPDATE  TT^F  RY  FLIGHT  TIME  OF  MATN  ROUND. 

TsTAFl T 

UPDATE  rPNTFR  OF  IMPACT  FOR  MAIN  ROUNDS. 

xrisRiASTGXR+XCAXBfXCI 

YCIsRPaSIGYRaYC^YH+YCI 

FAC  IS  A  FACTOR  TO  TNCLJDF  LAY  ERROR  TE  IT  IS  V- 
FACsI. 

enter  main  round  loop. 

.120 

9flS  Jr.J  +  1 

THin  TS  A  SWITCH  TO  TFLL  IF  THIS  TS  MAIN  ROIIMD  OR  MOT.  TH1021  IS  MAIN 
ROUND. 

IH1 0=1 
TH?0=0 

SAMPI  F  impact  point  for  THIS  ROUND 
TFfFAC.FO.O.^GOTOIOl 

IF  ROUND  HAS  HFF''*  SFNSFD  RFFORF,  EXClilDF  THE  LAV  ERROR. 

CALI  NRANS1 f Rl •RPI 
XPsSTGXLaRI 


Figure  7-1.  Source  Listing,  MAIN  Routine  (3  of  7) 
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YPsSTaYL*»? 

101  TAIL  MPAM^l (Rt •R?) 

C 

C  FTNr>  impact  POINT  OF  THTS  VATN  ROIINH. 

C 

XI  cXCI  fP1  *STC;xR  +  FAr*Y? 

Y1 SYC 1 ♦PPaSIRYR+FAC* YP 
C 

C  UPDATF  THF  TT'^F  FOR  THIS  MAIN  RIUNO, 

C 

CALL  NRAN^I (Tl , TP) 

TFf.I.FO.naOTniSR 

HWMsXMp*PXP(Tl*ST?) 

TF  f  I  F  .  AMT?  )  HiNM-AMTP 

TrTfH/yM^TFS 

onToi^n 

tS9  HWMsxMi  *FX«^  (  T  1  *ST  n 

TFfM'AJM*.!  F.  AMTn  HvrfMsAMTt 
TrT+HWM+TF 1 
C 

C  IS  TAPGFT  nFFTLAOF.n? 

c 

l#>n  TF  f  Y1  .IT.  YOFF)  GO  TO  1 
C 

C  IS  THIS  ROIIMO  WTTHTN  VUL  NF  R  AR  T  L  T  T  Y  MATRIX? 

C 

TFfXl^lT.AXf  1  DGnTOI 
TFf XI .GT. AX (MX+1 } IGOTOl 
TF  f  Y1  ^1  T,  AYd  )  IGHTni 
TF(Y1.GT.AY(NY  +  n)GnTni 

c 

c  vgHFRF  TG  THIS  ROlINn  TN  THF  Vt  JL  NF  R  A  R  T  L  T  T  V  MATRIX? 

C 

.11=1 

5  .Msllfl  , 

IF  f  Yl  .GT.AYf  jn  IGOTO-^ 

.?1  =.n  -1 

I?=1 

a  .fPrlP^I 

TF  (  X  1  .GT  .  AX  (.IP)  IGOTOa 
.IP  =  .IP-1 
C 

C  TS  THIS  ROHNO  A  HIT? 

C 

TFfTKf.II, , IP). LF.O  IGHTni 
C 

C  TF  THTS  TG  A  SPOTTING  RJFl.F  HIT  M^OATF  Mi|MRFR  OF  ROUNDG  TT  mOK  AMD 
C  START  FTRTMG'THF  matN  ROllNlD. 

C 

TFfTHIO^FQ.O)  TMM=THM+THM 
TF ( TH1 0 .F Q . n ) GOTnPna 
C 

C  TH  this  hit  a  RFLIARLF  Rnn\'n? 

c 

TF  f RAMHOMf  1  ) .GT .RFL )G0TnQ7 

r 


Figure  7-1.  Source  Listing,  MAIN  Routine  (4  of  7) 
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C  TF  TH^S  HTT  fS  NOT  THE  FIRST  HTT  GO  TO  5S 
C 

TFf  THI-^RT.OTGnTOSS 
C 

C  UPnflTE  THF  PRB&HTLTTY  ARRAY  ANO  TT^F  ARRAY  FOR  THIS  HTT. 

C 

THHsI 

.noz.i 

nnsi Ks.M  O.NROS 
51  7rK .^Ta7fK.6)>l . 

MOrP.  +  T/?. 

TFf  JlO.GT.ftnOOTnSS 

onisi.niaJi  n,f>i 

151  AKTlfJn.#>)aAKTL(I11.^)  +  1. 

c 

c  nOFS  THIS  HTT  HAVE  A  RFLTABI.F  FM7E? 

C 

55  TFfRANnnHn  ).GT.RFLFTr,nTq97 
THPOsl 
C 

C  TS  THIS  HIT  A  KTIL? 

r 

HWPzRANnnMd) 

KPPsTKT  J1  •.!?> 

MW1 =0. 

c 

C  9  OTFFFRFNT  KTIL  CRTTFRTA  ARF  mMSTDFRFD. 

C 

nOR^KPlzI.U 

TFcrrKpn.RT.o.TGnTnR^ 

TF  f PK (KPP.KPl ) .LT .MWPTnOTOR? 

C 

r  TF  THTS  HTT  TS  A  KTI.L  FOR  THIS  CRTTFRTOM  hpOATF  PRORARILITY  ARRAY  ANO 
’  C  TTMF  ARRAY.  * 

C 

.I10  =  P.  +  T/P. 

TF  f  ji  o.nT.hnnnTOi  sp 
oni  R-^.n  1  s,n  0.81 

15^  AKTI  f.fl  1  .KPnsAKTl  rjll  •KPI  )f1  . 

ISP  OnaPKPOs.f  ,  NROS 
9P  7rKPn,Kpn  =  7rKPo,Kpn>i - 

rfKPllrl. 

r,(un9^ 

RP  HWlst. 
q\  rOMTNIlF 

c 

C  TF  THTS  RnuNO  '^AS  NOT  A  KTIL  FOP  All  fPTTFRTA  GO  FTRF  ANOThFR  ROMNO, 

C 

TFfHi^M  .GT.O.inOTflPT 
TFfPASSM.r.T.o.  )r>nTnP7 

c 

C  TF  THTS  ROilHO  TS  A  KTI.L  FOP  All  rPTTFRTA  GO  TO  MFVT  SA^PL^. 

c 

GOTOqp 

c 

C  TF  THTS  TS  THF  v-iATN  RH'INH  r,n  Tn  pn  ^ 
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c 

1  TFf  TH1  o!,NF.0)f;nTO?0  5 

c 

C  TF  THF  SPOTTING  WTFLf  HAS  SHOT  MDWF  OR  fODAL  TO  THF  SPFCIFTFD  POUNOS. 

C 

TFfTHlI .GF.NROSt IGOTOPOh 
C 

C  y^AS  THTS  ROlINO  SFNSFD? 

C 

TFfRANnnwr 1 ),gt.pp  sn goto? on 
C 

C  rnRRFCT  n  tn  parallflogram  fashthn. 

c 

TALI.  NRAN^I  fRI  .RP) 
xrisxri«yMxc-Ri*STGxsi 
YCl =YC1-Vt ♦YC-RP*STGYS1 
FACrn. 

C 

c  FTRF  ANHTHFR  SPOTTING  ROUND. 

C 

GOTO^OO 

C 

C  CnRRFTT  THF  TFNTFR  OF  TARGFT  OtU  OF  CFNTFR  OF  IMPACT  BFCAUSF  IT  WILL 

C  BF  ADDFD  TN  WHFN  MAIN  ROUND  TS  FTRED. 

C 

?oa  xri=xri-yr 

Yfl SYCI -YC 
GOTOpni 

P0f»  THIPsIHiPtt 

IHI  ^=THM>THn 

c 

C  IF  ALL  OF  SPOTTING  RTFLF  RO  INDS  HAVF  BFFN  FTRFD  AND  ROUND  IS  NOT  WITH 
C  TN  IP  TNCMFS  OF  TARGFT  GO  TO  NFXT  SAMPlF  OTHFRwTSF  FIRF  MATN  ROUND. 

“  C  .  • 

TFfAXfl  I-l?..GT.xn  GOTORR 
|FfAyfNyif1?..l  T.X!  )  GOTORR 
TF(AYrNYWt?..LT.Yl  )  GOTORR 
TFCAYfl  l-lP..GT.Yn  GOTORR 
GOTOpn/J 
C 

c  TF  THF  MTSS  OF  THF  MAIN  ROijMD  TS  NOT  SF'NSFD  FTRF  ANOTHER  ROUND  WITH 
C  NO  rORRFCTTON. 

c  * 

PC'S  TF  (RANDO^^n  )  .GT  .  PROnSI  GO  TO  R7 
C 

C  CORRFCT  ri  FOR  NF.XT  ROUND  USING  P  AR  AL  L  FL  OGR  AM  METHOD. 

C 

CAL  L  ^*RAN  f  RI  ,  RP  ) 

xrisxrt-xi+xr-Ri *stgxs 

YCi  =Yri -Y1  ♦YC'*RP*STGYS 
FACrO. 

R7  rONTTMl'F 

TFfPASSN.!,F.O.IGOTORB 

C 

C  CHFCK  FU7F  ACTUATION, 

C 
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TFf TH?0.FO.0)GnTO?0? 

C 

C  PFRSnNNFL  rASUALTirS  COMPOTFO. 

c 

SnvirAS  sSDMCAS  ♦(PK(K??,S)/PASS'M)*(PASSN*Sll'^CAS  ) 

\(\?  7(.t,S)s7(J.S)+SUMrAS 
Qft  rONTTNUF 

c 

C  FTPF  NFXT  main  POUND. 

C 

TF( J.L  T.NPDS^GOTOPflS 

C 

C  Gn  TO  NFXT  SA^^PLF• 

C 

99  rONTTNUF 
« 

c 

C  SAMPl  TNG  TS  CO'^PLFTF. 

C 

wPTTF(fe,7) 

7  FOP'^ATf //tPW  NPDS  PK/) 

C 

C  CO'^'PUTF  MFAN  PK  AND  PH  TAPGFT  DATA  AND  /*PTTF  PFSULTS, 

C  , 

SA^P-nSAMP 
SAMPTsI  ,/SA'^P 
DO  9^  T?=1 .NPDS 
DOPI  T  =  1 

91  7f TP.Tl=7fT?.T)*SAMPT 
«7  wPT  TF  f  7  (  T?.  T  )  .  Tsi  .8) 

9^  rnNTTNIlF 

THt  ^sTHI ^/NSAMP 

c 

C  TH1P  TS  NHMMFP  DF  TT'^FS  FTPFD  ALL  SPf3TT|MG  ROUNDS.  IHM  IS  AVFPAGF 
‘C  NUMRFP  OF  SPOTTING  POUNDS  FTPFD. 

c 

WPTTFf^.PSn  THIP.THM 
PS1  FOPMATTPISI 
C 

C  rO^^PUTF  mfan  PK  AND  PH  TT'^F  DATA  AND  ;a/PTTF  PFSUlTS. 

C 

WPTTFf8.#»n 
DOI SaT=1 .^1 
DOiSajsi.f* 

ISa  AKTl  f  I  ..TIsAKTLf  T  .  II*SAMPT 
DOISS  T  =  1  .S 

'nRT  TF  f  1  SD  W  Ak  T  L  n  .  I )  .  T  =  1  .  6  t  ) 

ISf,  FORv^AT  f  1  OFft.SI 
WPTTFT^.1S7) 

1S7  FnP^'ATflMn) 

1SS  COM  TIN' IF 

c 

C  rOMPHTF  AVFPAGF  NO.  OF  Pns,  FTPFD. 

C 

TFf  NOPT1  .F^.'^IGOTOP 
TAIL  KTI.PAT 

r 

C  GO  TO  PFAD  NFXT  DATA  SFT. 

C 

GOTO? 

6  FOPMATf TS.  9Fl%. 71 
hi  FORMATTIHI) 

h?  FnPMAT(//P1H  NUMRFP  OF  SAMPLFS  =  ,Tfc./tHn 
FND 

Figure  7-1.  Source  Listing,  MAIN  Routine  (7  of  7) 
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FIIMCTTON  CMF(X) 

C**********« 

c  PURPOSF 

,C  C^-l^PllTFf^  ruvtULATTVF  hi)Q^&L  DTSTRTBUTTIN  VALUFS  FOR  CFLL  PH  tiSFO 
'c  TM  rinSFO  FORM  snuJTTOM. 


F  =  0. 

AX=AnSfX)*.70710hR 
TFC  AX.r^F.S.  )aOT1^ 
Fs.S-.S*FRF( AX) 

TFf X.np.O.IFsl .-F 

rNF  =  F 

RFTIJRN 

FMD 


Figure  7-2.  Source  Listing,  Function  CNF(X) 
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r»r>o  r»r>n  nr»n  n  n  n  r>r»r»  noo 


.SURt^DMTTNF  T^^PUT 


«****«*« 


C  PURPOSF. 

C  RFAO  AMR  '.VRTTF  INPUT  DMA#  CONVFRT  UNITS  OF  INPUT  DATA,  ANO  vSFT  UP 
THF  TAPnFT  RPin  ANO  PK  ARRATS, 


niM^NSTON  PKf^OOO  •SI,  AXtPni  ),AY(S1  )  •CASF(iiO)fZftP#9)#7H6n 

rOMMHM  PK.AX.AY,SXR,Sy».SXR,SYH,SXS,SYS»PRS»Xfl#YR,XC#YC,RFL#TFl, 

J  TFS, XM1 , XMp, YnFF.NOPTI , MRDS . NX . NY , FL T , S T t •ST?,AMT1,AVT?,IK(S0#?00) 
?*7, 71  •  AKTL  (81  *9)  .SIRXL,  SiaYL#RFI.F 

^•S  XS1.S  YS1.S  XR1,S  YRt,S  XR1.S  YR 1 • PRS 1 • XR t • YB 1 • FL T1 • TF 11 » 
aTFSl  •  XMll  ,  XM?t  •  AMTI  1  ,  AMIPI  ,  ST  11  •ST.?l  •  SIGXLl  .SIGYL  1  •NRDSl  .PASSN 

RFAD  AND  WRTTF  TITLF  CARO. 

PFAnfS,?n.FN0s6l  CASF 
?ft  FORVATfaOAPI 

WRITFfh.lSI  CASF 

lA  FORMAT  (  IHl  •ilOAP) 

RFAn  main  round  DATA. 

RFADfS, 1 •FND=8)N0PT 1 • NROS • WOP T • XT • YC • SXS • SYS . SXR • SYR • SXR • SYR . PRS # R 
IFl  . YDFF.XR, YR.FLT, TF1 • TFS*  XMl . XM?, AMT t • AMT? , ST  1 • ST?, SIGXL • STGYL # 
?RFIF,PASSN 

CHFCK  FLAG  TO  RFAO  SPITTING  ROUND  DATA. 

IF (NOPTllNF.PIGOTnPO 1 

RFADfS»?na)  NRDSl  .SXS1#SYSl.SXRl,SYHl,SXRlfSYRt0pR5>l,XRl»YRlfFLTK 
tTF11,TFS1,XMM#XM?t,AMTn»AMT?1pSTn#STPl.STr.XLt»STGYLl 
?n«  FORMATfTS,SX,7F10.?./flFin.?,/RP10.P) 

NUMRFR  of  main  ROUMDS  FTRFD  PFR  SAMPLF  SPFCTFTFn. 

?01  TFfNRDS.L T. 1 0)  NRDSslO 
TFfNRDS.GT.lRI  NRDSslO 
TFf NOPT1 .FQ. NRDSsl 

WRTTF  MAIN  ROUND  INPUT  DATA. 

WRTTF r^,?l NOP T1 /NR^S, NOP T 

?  FORMAT!  ////////7//P1H  ST^'UI.ATTOM  s  .TS^/POU  NUMRFR  OF  R 

10UNDS  r  .TS*/flH  NOPT  =  ,TS) 

WRITFf^.n 

*5  FORMAT! //?RX#?aH  HORT7nNTAL  VFRTTCAL  ,//SlX,7M  M^fPRS,  SX,7H  ME 
ITFRS,/) 

WOTTF  ff>.q)SXS.SYS,SXM,SYR,SXR,SVR,XR, YR,XC» Yf •STGVL.STGYL 
a  FORVAT!?#>H  SENSING  FRROPS  !STD  nrV)  •  ?F  t  ? .  y  , /?41H  VARTARLF  BIAS  (ST 

ID  OF VI  •  px •  ?n ?. • /pan  random  frrors  (sto  df vi  •  px , pfi  p. /i ,  / 1 1 h  fixfd 

P  HI  AS,  1  RX  •  PF1  P,  41  ,  /  1  1  M  ATM  rFNTFR  ,  1  SX  ,  PF  1  p  .  41/ ,  1  ON  lAY  FRR0R,1^X, 
SPF1P.41/I 

WRTTF !S* 7 )PRS,RFL, YnFF,RFLF 

7  FORMAT!  iQM  PRORARTLITY  OF  SFNSTNG  MISSING  RDU^fD  r  ,  F  1  ?  .  6  , /?  RFLIA 
1BTI.TTY  OF  ROUND  =  VFRTTfAi.  DFFTI  ADFD  ^O^RnT^fATFS  r  , 
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?F1?,8.7H  TNCHFS/.P^H  RFL.  T  AR  T I.  T  T  Y  OF  FU7F  = 
wmTFff^^l^ilFLT.TFI  .TFS,X^^!  .¥MP,ANflT1  .aMT?,ST1  .ST? 
q  FFIRWATf 18H  TT'^F  0  A  T  A  f  SF  C  S )./// 1  RH  TTMF  OF  FLIGHT  =  .FI?. 6 
1.///?RX.inH  tST  ROIIMO, 9K , 1 RH  SURSFOUFMT  RnS..//?SH  FTXFO  MFHTAN  TI 
?MF  .FTP.6. 10X,F1?,R./?SH  VARIARLF  VlFpIAN'  TIMf  .F1?.R.10X. 

■^F1?,R./?SH  MTMIMOM  TTMF  ,  F  1  ?  .  6 ,  t  0  X  ,  F  1  ?  .  R . /?SH  SLOPE  OF 

an T S TP T RUT  TON  ,F1?.R. 1 OX.FlP.Rl 

CONVFPT  MATM  PO.  OFl  TVFPy  DATA  FROM  ^^FTRJr  TO  FNGLTSH  STA^f^ARD  UNITS. 

fiXSrSX.S*XR.^7 
SYS=SYS*XR.X7 
SXR=SXR*XR. X7 
SYR  =  SYR*  ^9. X7 
SXR2SXR*XQ. X7 
SYRrSYR*X9.X7 
STGXL  sSTGVL 
STGYl  =STGYL<>  ^9.'i7 
XR  sXH 

VR  SYR  <k  X9  .  \7 

xc  =  xr  ^ 7 
YC=Yr*X9.X7 

CHFOK  FLAG  FOP  WRITING  SPOTTING  PTFLF  ROUND  DATA. 

TF  f NOPTl .NF .?)GOTn?n? 

RFLF1 St  . 

RFLIst  . 
wRTTFTR.IR) 

WRTTFfR.POXl 

?OX  format ( ///?nH  SPOTTING  PTFLF  DATA.//) 
rtPTTFTR.X) 

wPTTFTR.alSXSl.SYSl ,SVR1 .SYR1 .SXRt.SYRl .XRl.YRlfXC.YC.SlGXLl.STGYL 
11  • 
WPTTFfR.7)PPS1 .PFL1 . YOFF.PFLFt 

WPT  TF (R. 1 PIFLT  t . TFll .TFS1  . XMl  1 , XM?t . AMTM . AMT?1 ,STt 1 .ST?1 

rONVFPT  SPOT  pr>.  DFLTYFRY  data  FRO^  MFTPTC  to  ENGLISH  STANOAPn  UNITS, 

SXS1 sSXS 1  *  XR  .  ^  7 
SYSl sSYSt * ^P. W 
SXRIsSxni^XP.'^? 

SYR  1 sSYRl  A X9. X 7 
SXRt  sSXRI  *X9.  ^7 
SYR1 sSYPI *  ^9. X7 
STGXI  1  sSTGyi.  1  *  S9.  W 

STGYL 1 sSTGYL 1  *^9.^7 

XR 1 SXRt  A  ^9 . X7 
YRI sYRl A  XP.  W 

CPFW  PMOFTrTFNGY  TT^F  FACTORS  SFT. 

TFfSTl t .FQ.n.I  Sit ls.a9R^ 

TF(ST?1  .FO.O.  )  ST?ts.a9RX 
?ft?  TF  (ST  1  .FO.O.  )  STls.fjQSS 
TFTSTP.FQ.n.I  ST?s./]9S^ 

10?  TF(X.GT.999, )  RFTURN 

C  IMVFRT  TAPGFT  COOROTMATFS  INTO  TO  K  T  M.  ARRAY  SI.IRSCRTPIS, 
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c 

C  CHFfK  INPUT  C()NT»?OL  FLAG  FOR  USING  PRFVIOns  DATA*  SFT  UP* 

C 

IFfNOPT.KO.OIRFTUR^J 

C 

C  RFAH  TARGFT  DATA. 

C 

RFAHfS,  nWX.MV 
NYtsNY+l 

I  FnWMATf ^IS,SX.6FtO,?*/RFlO.P# /flFl 0.?, /SFI 0.?) 

IS  RFAn(S,St iXMTN.YMTN.OX.nV 

c 

C  SFT  VULNFRARTLTTY  CARU  LTwIT. 

c 

NC  =  '^n  00 

c 

C  SFT  UP  VIH.NFRARTL t lY  MATRIX. 

C 

AXf I IsXMTN 

nn  iR  TsP.NXi 

1 h  AX  f  T 1 =AX  f I -1  I ♦DX 
A  Y  (  n  =  Y  M I M 

nn  17  isP.MYi 

17  AYfnrAYf  !•!  l^DY 

C 

C  1NTTTALT7F  ID  KILL  ARRAY. 

C 

S8  nnsaisi.Nx’ 
nnsajri .ny 
Sa  TKfI.T)=0 
C 

C  FNTFR  VUI NFRARTLTTY  CAROLOnP. 

c 

nns?T=i ,ur 
MtlrS  , 

TFrPASSN'.I.F.O.)Nri=a 

c 

c  rHFGK  FLAG  FOR  STANOARD  ARRAY  SFT  UP. 

C 

TFfNOPT.FQ.I  IGOTO^O'^ 

C 

C  RFVFRSF  TARGFT  HRIFNIATinN. 

C 

RFAOfS.  >  X,  Y,  f  PK(  T  .  n  #.J=1  -  NCI  ) 

y  =  -x 

GOTn^np 

50*^  RF  Anf  s. -^os)  X,  Y,  (PK  ( T ,  Ji  •  »  =  i .  Mrn 

■^n/j  format  ( PF  7 .P,  , 'jF  n 

5ns  fhrmakpft.p.  A»X.SF7.n 

c 

C  rUFTK  FOR  I  AST  V '  )l.  NF  R  A  R  T  L  T  T  Y  CARn. 

c 

J1z(fX-XMTN)/nX)Fl  . 

JPs ( f  Y-YMT  N) /OY  T ♦!  . 

SP  TKf.JP,Jt)zT 
10  FORMATflMIT 
RFTURN 

Si  FORMAT CaFI 0 . 1 . IS) 
h  CALI  FXTT 
Fwn 


Figure  7-3.  Source  Listing,  Subroutine  INPUT  (3  of  3) 
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noo  ono  nr>r»  r>r»r>  omri  orinnon 


SDHROIJTTNF  KTLPAT 


PflRPdSF 

CiMPtITFS  AMD  WPTTFS  THF  AVFPAGF  NtlMf^FR  OF  MAIN  ROUNDS  FTRFD  PFR  CASE. 


********** 


***************** 


OTMENSTON  AKTLf  ^I.R^FC  11  •  r,  (  6n,SMKCiq  ,R),AA(  h  1  1  ,  PK  f 0  0  •  S 1  , 

1  Axfpol  1 .  AY  f  sn 

COMMON  PK.AX, AY.STGXR.STGYR.STOXR,  S T G Y B , S T GXS , S T G YS , S I G  ,XR,Y8, 

1  xr  ,  YC.RFl  .  TF  1  ,  TFS,  XM1  ,  XM?  ,  yOEF  ,  MOP  T  1  ,NR  •  NIX  ,  MY  .  FL  T  .  ST  1  ,  ST? ,  AMT  1 , 
?AMT?, TK (S0*?00) ,SMK,F , AK  TL , STGXL  #  SIGYL . WFLF 

^•STGXSI .STGYS1 *STGXB1,STGYR1 ,STGXR1 ,STGYR1 ,PRS1 ,X81 , YB1 ,Fl  Tt  *TF11 , 
iiTFSl  .  XMt  1  .XM?i  ,  AMT1  1  ,  AMT?1  .  ST  11  ,  ST?1  ^STGXL  1  ,  STGYL  1  •NROS1  .PASSN 


FNTFR  case  I  OOP. 


DO  ?  T 1=1  .6 
SMJKsSMKd  ,111 
AUL=?. 

AMRsNR 
.?  =  NR 


FMTFR  ROUND  LOOP  AND  COMPUTE  AVG.  RDS.  FIRED. 


DO?ST=?.J 

SMTK2SMTK*( AUL* (SMK ( T , T 1  1-SMK ( T-t , 1 1 1 1 1 
?S  AUL  =  AUL*1. 

DTD  ALL  SAMPlFS  REGISTER  A  KILL? 

TF  f SMK  CNR, T 1 1 ,GT. .RRRRl  GOTORR 

A«  =  S.*(SMK U-S, I  1 1 ♦SMK  f J-a, I  1  1 ♦SmkC T 1 1*SMK ( J-?, 1 1 1 4SMK ( J-1 , T  t 1 ) 

11 

XXsSMKf.I^  ,  1 1  1-Smk  ( J-S,  1 1  1 
^  TFrXX.LF.O. IGOTORR 

A7=AR/XX 

COMPUTF  AVG  ROS  FIRFD  TMCt.UDTMG  NON-K  T  t.  L  T  MG  SAMplFS  AND  EXTRAPOLATION  FACTOR 

AVs(A74ANRl*n  .-SMK  (.1  .TIIWSMTK 
GOTORfl 

WRTTF  AVG  ROUNDS  FIRED  PER  CASE. 

99  AVrSMTK 

RR  WRTTFr4».?Rl  AV,I1 

?h  FORMATfl^H  AVG.  MO.  RDS  =  ,F7.?.HH  CASE  =  ,TS,//1 
?  rONTTNUF 
RETURN 
END 


Figure  7-4.  Source  Listing,  Subroutine  KILRAT 
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iT  X 


SURMniJTT»^r  NRAN^I  (X!  ,  Y?1 


C  PURPHSF 

C  r.FNFRATFS  RANDOM  NHR'^AL  HFVTATFS. 

C 

C  FUNCTTON  RAMOn'^  RFOIJIRFD. 

C***************<>****^*****^********^********^******'*****^*****^**************** 

C 

Y^5SQRT(-?.*ALOG(RANnnw(l ) ) ) 

Xy  =  S,PA-^t  RS^0  7?*RA.Nn0M(  !  ) 

XP  =  X  ^♦STNf XO) 

y 1 =x  ^*rns(xai 

RF  TURN 
FNO 


Figure  7-5.  Source  Listing,  Subroutine  NRAN31  (XI,  X2) 


FIINCTTOM  RAMOHMCN) 

C*******4*A******************«A***A**4*ft***********«******A**'**A**A«*A*******«*4 

puRpnsF 

r^FNFRATFS  UNIFORM  RANDOM  NUM^pj^s  RFH»’FFN  0.  AND  t. 


c 

rOMMON/  RANDM/TX 
TXrlX*f>SS\R 
TFf  IX)S,^.(S 
TX  =  TX  +  PW4';U8^6a74'l 
RUM  rXX 

t  SUM  =  SUM*'.aftS681  "^F-R 

R AMDOMsSUM 
RF  TURN 
FNO 


Figure  7-6.  Source  Listing,  Function  RANDOM  (N) 
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r»or»  nor>  nr>o  non  ooo  nno  ooo 


- SIJRROUTTNF  SSMOT  ( PC  H ,  PC  W  H 

C 

C  PMRPnSF 

C  Cl'^PHTFS  pH  TAWiFT  FOR  A  STNRLF  '^ATM  ROUND  AND  FOR  A  STNGLF  wa|W  ROUND 

C  GTVFN  A  SFNSFD  MISS  RY  USING  A  CinSFO  FORM  SOLUTION. 

C 

C  FUNCTION  CNF  RFOUIRFD, 

C*i»A**A**»A****^>***************************************************************^ 

c 

;  nTMFNSTON  PKf^OOO  . 8  I . A  X ( PO IT , A Y ( S 1  I , ? C 1 P , P ) . 7 1  f 6 1 I 

COMMON  PK,AX,AY,SIGXR,SIGYR,SIGXR,  S 1  GYP • S I GX S , S I GYS • PRORS * XB , Yfl , 
tXC,  YC^RFL.TFl.TFS.XMt  ,  XM;>.  yOFK  ,  NQPT  1  .  NROS  ,  NX  .  N  Y ,  FL  T  •  ST  1  ,  ST? ,  AMT  !  , 
?AMT?.TKfS0  .?0n),/,7l ,AKTI ( h t , 0 1 • S T G X L • S I G Yl , RFL F 

STGXS1  ,  SIGYS1  ,  STGXR1  ,  SIGYHI  ,  STGXR1  ,  SIGYRI  ,PRS1  ,  XR  1  ,  YR  1  ,  FL  T  t  •  TF  1  1  , 
UTFS1.Xwn.yM?i,AMTtt,AMT?l,STtl,ST?l.SIGXLl , SIGYL 1 ,NRnSl .PASSN 
PCH  =0. 

PCH1 20. 

COMPUTF  impact  sigmas  FIR  A  STMGLF  maTM  RTUMD. 

SIGXrSORTf  STGXR**?4'SIGXR**?^STGXL**?) 

SIGYrSQRT  f SIGYR**?  +  STGYR**?>STGYL**?I 
STGArSDRT  f STGXS**?^?.*STGXR*  *?) 

STGRrSORT  (SIGYS**?  +  ?,i*STGYR*  *?) 

CHFCK  FOR  SPOTTING  ROUND  DATA, 

IFCNOPTI .NF.?)GOTO^ 

MODIFY  impact  sigmas. 

STGXzSORT  f  SIGX<^*?f  (  ,  S*  f- S  I GX  S 1  ♦  S  I G  XL  11  )  *  *  ?♦  S  I  GXR 1  *  *  ?  I 
SIGYrSQRT  CSIGY*  (  .S*  (SIGVS  t  ♦STGYLn)**?  +  SIGYPl**?l 

FNTFR  LOOP  TO  SFLFCT  CFLL  VFRTTCAL  COORDINATFS. 

3  DOPJrl.NY 

Y?2AY(J  +  n-YC 
YlrAVf.f  )-YC 

CHFCK  lARGFT  DFFILADFD  rnNOITIOMS. 

TFT  AYM^n.LF.YDFn  GOTO? 

■  TFfAYfJ  l.LT.YDFFl  YIzYDFF-Yr 

FNTFR  1  OOP  TO  SFlFd  CFLL  hoRT/Omtal  COORDINATFS. 

001 l2l .NX 

CHFCK  THAT  CFLL  IN  TARGFT. 

TFflKf.I.n.LF.O  IGOini 
XlrAXd  I-XC 
X?2AXf If  1  )-XC 

c 
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C  CO^PUTF  ANH  CFLU  ^H*S. 

C 

PA  sCtNPf  f  )(P-yB)  /STGXl  -CNF  (  f  X  t  - VH  ) /S  T  G X  )  )  *  ( CNF  C  (  Y?-YR)  /STGY  ) -CNF  f  ( 
1  Y1-YH)/3TGY1 ) 

PA1  rfCNF  f  X?/S-TGA)-CNF  f  X  1  /  STGA  )  )  *  (  CNF  (  Y<?/S  I GR  > -CNF  (  Y!  /  STGB)  ) 

PCM  sPfH  fPA 
PCHJrPCHI tPAI 
1  rnNTTNtlF 

?  CDNTTNIIF 

PFTURN 
FNO 


Figure  7-7.  Source  Listing,  Subroutine  SSHOT  (PCH,  PCHl )  (Page  2  of  2) 
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APPENDIX  A 


Appendix  A  contains  the  entire  vulnerability  card  input  deck 
the  sample  case  presented  in  Section  VI. 


A-1 


?ti .  0 

flo  .  0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0,000 

.0 

00.0 

0.000 

O.OOO 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

ao ,  f) 

^p.o 

0.000 

0.000 

0.000 

0.000 

0.000 

o.ooo 

0.000 

0.000 

■\6.o 

^P.O 

0.000 

0.000 

0.000 

o.ono 

0.000 

0,000 

0.000 

0.000 

^P.O 

.987 

.998 

.998 

,  798 

o.ooo 

0.000 

0.000 

O.OOO 

^fl.O 

iP.O 

o.ono 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

;?y  .0 

^?.0 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.0 

^P.O 

0.000 

0.000 

0.000 

0.000 

o.ooo 

0.000 

0.000 

o.ono 

o.n 

^P.O 

0.000 

o.ooo 

o.ooo 

0.000 

0.000 

o.ono 

0.003 

0,00  0 

-0  a .  0 

^P.  0 

0.000 

0.000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

-00,0 

^P.O 

0.000 

0.000 

0.000 

o.ono 

0.000 

0.000 

0,000 

o.ooo 

-s?.o 

^P.O 

0.000 

0.000 

0.0  00 

0.000 

o.ooo 

0.000 

0.000 

0.000 

00 . 0 

PA.O 

.987 

.998 

.998 

.  798 

0,000 

0.000 

0,000 

0.000 

00.0 

Pfl.o 

.9fl  7 

.998 

.998 

.  798 

0.000 

o.ooo 

o.ooo 

0.000 

X^>.0 

PR.n 

.987 

.998 

.908 

.798 

o.ooo 

0.000 

0.000 

o.ono 

^;>,o 

PH.O 

.9A7 

.  998 

.998 

.798 

0.000 

0.000 

0.000 

o.ooo 

?R.O 

PA  .  0 

.987 

.998 

.908 

.798 

o.ooo 

0.000 

o.ooo 

0.000 

?O.0 

PA.O 

.98  7 

.998 

.998 

.798 

0.000 

0.000 

0.000 

0.000 

?o.o 

PA.O 

.987 

.998 

.998 

.798 

0,000 

0.000 

o.ono 

0.000 

1R.0 

PA.O 

O.ono 

0.000 

o.ono 

0.000 

0.000 

o.ooo 

0.000 

0.000 

1  ?.o 

PA^O 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

R^O 

PA.O 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000  ’ 

0,000 

0.0 

PA.O 

0.000 

o.ono 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0.0 

PA.O 

o.ooo 

o.ooo 

0.000 

0,000 

0.000 

0.000 

0.000 

0.000 

-0.0 

PA.O 

0.000 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

o.ooo 

-fl.O 

PA.O 

0,000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-Ofl.O 

PA.O 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0.000 

0.000 

-s?^o 

PA.O 

0 .000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

00.0 

PO.O 

0 .000 

0.000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

O.OOQ 

00.0 

po ,  0 

.98  7 

.998 

.998 

.798 

o.ooo 

0.000 

0.000 

0,000 

00.0 

PO.O 

.987 

.998 

,Q98 

.798 

.987 

.998 

.998 

.798 

po  ,  0 

.987 

.99  8 

.998 

.798 

.987 

.998 

.998 

,798 

0 

Po  .  0 

0.000 

0 .000 

0.000 

O.OOO 

1  .000 

1.000 

l.OOO 

l.OOO 

PO.O 

1  .000 

1  .000 

1  .000 

1  .000 

1.000 

l.OOO 

1  .000 

1.000 

?n .  0 

PO.O 

1  ,000 

1  .000 

1.000 

t  .000 

1  .000 

1.000 

l.OOO 

1,000 

?o^o 

PO.O 

0.000 

0.000 

0.000 

o.ooo 

1  .000 

1  .000 

l.OOO 

l.OOO 

1  8,0 

?0 . 0 

O.OOO 

0,000 

o.ono 

0.000 

0.000 

0.000 

0.000 

o.ono 

1P.0 

PO.O 

O.ooO 

o.ono 

0.000 

0.000 

.987 

.998 

.998 

,798 

fl.O 

PO.O 

0.000 

0.00  0 

0.000 

0 . 0  0  ii 

0.000 

o.ooo 

o.ooo 

o.ooo 

o.n 

PO  .  0 

O.ooo 

0.000 

0.000 

0.000 

0.000 

0 .000 

0.000 

o.ono 

0.0 

PO.O 

0.000 

0 .000 

o.ooo 

o.ooo 

0.000 

0 .000 

o.ooo 

0.000 

-0.0 

PO.O 

0 .000 

0.000 

0.000 

o.ono 

0.000 

0.000 

0.000 

0.000 

-fl.O 

PO.O 

0.000 

o.ooo 

0.000 

0.000 

0 .000 

0.000 

o.ooo 

0.000 

-i?.n 

po  .  0 

o.ooo 

0 .000 

0.000 

0.000 

0,000 

0.000 

0.000 

0  .  noo 

-18.0 

PO  .  0 

O.ooo 

0,000 

0 .000 

0  .  OOO 

0.000 

o.ooo 

0.000 

o.ono 

-PO.O 

PO.O 

o.ono 

0 .000 

o.ooo 

O.ooo 

0 .000 

o.ono 

o.ooo 

0  .  n  0  0 

-■\P.  0 

PO.O 

0.000 

o.ooo 

0 .000 

0.000 

0.000 

.  0.000 

o.ooo 

0.000 

-*^8.0 

?0 . 0 

0 .00  0 

o.ooo 

0.000 

0.000 

,987 

.99  8 

.99  8 

.  798 

•  0  '4 . 0 

PO .  n 

0 .000 

0.000 

o.ooo 

0 , 0  00 

0 .000 

0.000 

0.000 

0 .  non 

Ofl^o 

PO.O 

0.000 

0.000 

0.000 

0.000 

0.000 

0,000 

0,000 

0,000 

00^0 

PO.O 

1 .000 

1  .000 

t  .000 

1  .000 

0.000 

0.000 

0,000 

0,000 

00^0 

PO.O 

0.000 

0.000 

0.000 

0.000 

1,000 

1,000 

1,000 

1 .000 

^8.0 

PO.O 

0.000 

o.ono 

0.000 

0.000 

1 .000 

1 .000 

1  .000 

1  ,000 

PO.O 

1 .000 

t  .000 

1  .000 

1  .000 

1 .000 

1,000 

l.OOO 

1.000 

?fl,0 

PO.O 

0.000 

0,000 

o.ooo 

0.000 

1.000 

1.000 

l.OOO 

l.OOO 

?o^o 

PO.O 

1 .000 

1  .000 

I  .000 

1  .000 

1  .000 

1  .000 

1 ,000 

l.OOO 

PO.O 

PO.O 

0.000 

0.000 

0.000 

0.000 

t  .000 

1 ,000 

1.000 

1  .000 

1  8.0 

PO.O 

.987 

.998 

.998 

.798 

.987 

.998 

,998 

.798 

IP^O 

PO.O 

.987 

.998 

.998 

.798 

.987 

,998 

,998 

.  798 

fl.O 

PO.O 

0.000 

0.000 

0.000 

0.000 

.987 

.  .998 

,998 

.798 

0.0 

PO.O 

.987 

.998 

.998 

.798 

.987 

.998 

.998 

.790 

0.0 

PO.O 

0.000 

0 .000 

0.000 

0.000 

.987 

,998 

,998 

.798 
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n 

?n,o 

0.000 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0,000 

?o.n 

o.ooo 

0.000 

0.000 

0.000 

,987 

,998 

.998 

.798 

-i?.n 

PO.O 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

O.OOO 

-t6^0 

;>n.n 

0.000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-?0.0 

?0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

?0.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

?o.n 

0.000 

0.000 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

1 

x» 

• 

o 

PO.O 

.9ft7 

.998 

.998 

.798 

0.000 

0.000 

0.000 

o.ooo 

•■^6^0 

po.n 

.9A7 

.998 

.998 

.798 

0.000 

0.000 

0.000 

o.ooo 

-uo.o 

PO.O 

0.000 

0^000 

0.000 

0.000 

.987 

.998 

.998 

.  798 

-aa^o 

PO.O 

0.000 

0.000 

0.000 

0.000 

.987 

.998 

,996 

.  798 

s?^o 

18.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

aA^o 

18^0 

.9A7 

.998 

.998 

.798 

0.000 

O.OOO 

0.000 

0.000 

aa.o 

18.0 

0.000 

0.000 

0.000 

0.000 

.987 

.998 

.998 

.798 

ao,o 

18.0 

1  .000 

1 .000 

1  .000 

1  .000 

1  .000 

1 .000 

l.OOO 

1  .000 

18.0 

1 .000 

1  .000 

1  .000 

1.000 

1  .000 

1,000 

1  .000 

1  ,000 

18^0 

1  .000 

1  ,000 

1.000 

1  .000 

1  .000 

1  .000 

1  .000 

1  .000 

?fl.O 

18.0 

1 .000 

1  .000 

1  .000 

1 .000 

1  ,000 

1 .000 

1  .000 

1.000 

?a^o 

18.0 

.987 

.998 

.998 

.798 

.987 

,998 

.998 

,  798 

PO.O 

18.0 

.987 

.998 

.998 

.798 

.987 

.998 

.998 

.798 

1h,0 

18.0 

.987 

.998 

-  .998 

.  798 

.987 

.998 

.998 

;798 

1  ?.o 

18.0 

1 .000 

1  .000 

1  .000 

1  .000 

1 .000 

1,000 

1,000 

1.000 

R.O 

18.0 

.987 

.998 

.998 

.798 

1.000 

1.000 

1  ,000 

1.000 

a.O 

18.0 

.98  7 

.998 

.998 

.798 

.987 

,998 

.998 

.798 

0.0 

18.0 

,987 

.998 

.998 

.798 

.987 

.998 

.998 

.  798 

-a.o 

18.0 

.987 

.998 

.998 

.  798 

1  .000 

1 .000 

1  .000 

l.OOO 

-fl.O 

18.0 

.  .987 

.998 

.998 

.798 

1 .000 

1  .000 

1  .000 

1  ,000 

-t?.o 

18.0 

.987 

.998 

.998 

.  798 

.987 

.998 

,998 

,798 

•18.0 

18.0 

.987 

.998 

.998 

.798 

.987 

.998 

.998 

,798 

-?0.0 

18.0 

.987 

1  .000 

1.000 

.798 

0.000 

0,000 

0.000 

0,000 

•?a  .0 

1  8^0 

.987 

1  .000 

1  .000 

.798 

0.000 

0,000 

0.000 

0.000 

18.0 

.987 

1  .000 

1.000 

.798 

0.000 

0.000 

0.000 

o.ooo 

18.0 

0.000 

o.ooo 

0.000 

0.000 

o.ooo 

0.000 

0.000 

0.000 

•^6.0 

18.0 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

n  .000 

•ao.o 

18.0 

.987 

1  .000 

!  ,000 

.798 

0.000 

0.000 

0,000 

0.000 

•aa.o 

18.0 

0.000 

0.000 

0.000 

0.000 

.987 

.998 

.998 

.798 

-afl.o 

18.0 

0.000 

1 .000 

1  .000 

0.000 

.987 

,998 

.99  8 

.798 

-s?.o 

18.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-78.0 

18.0 

0.000 

1  .000 

1  .000 

0.000 

0.000 

1  .000 

1  .000 

0,000 

-AO.O 

18p0 

0.000 

1  .000 

1  .000 

0.000 

0.000 

1,000 

1.000 

0.000 

SP.O 

IP.O 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

aR.o 

1P.0 

0.000 

0.000 

0.000 

0.000 

1  .000 

1 ,000 

1  .000 

1.000 

aa.o 

IP.O 

0.000 

0.000 

0.000 

0.000 

1 ,000 

1  .000 

1  .000 

1  .000 

ao.o 

IP.O 

.987 

.998 

.998 

.798 

1.000 

1  .000 

1  .000 

1 .000 

u.o 

IP^O 

.987 

.998 

.998 

.  798 

,987 

,998 

.998 

.798 

^?.o 

IP.O 

.987 

.908 

,9t^8 

.798 

1  .000 

1,000 

1 .000 

1  ,000 

?A.O 

1?;.0 

.987 

.99T 

.998 

.791 

,987 

.998 

,990 

,798 

?a.o 

1P.C 

.98  7 

.  9"^(i 

.998 

.  798 

.987 

.998 

.998 

.798 

PO.O 

IP.O 

.987 

.998 

.998 

.798 

1  .000 

1 .000 

1  .000 

l.OOO 

18.0 

IP.O 

.987 

,998 

.998 

.798 

.987 

.998 

.998 

.798 

1?.0 

IP.O 

.987 

.998 

.998 

.798 

1 .000 

1,000 

1  .000 

1 .000 

fl.O 

IP.O 

1  .000 

1  .000 

1  .one 

1  .000 

1  .000 

1 .000 

1  .000 

1  ,000 

a^o 

IP.O 

.987 

.998 

.998 

.798 

1 .000 

l.OOO 

1 .000 

1  .000 

0.0 

IP.O 

.987 

.998 

.998 

.798 

1  .000 

1.000 

l.OOO 

1.000 

-a^o* 

IP.O 

.987 

.998 

.998 

.798 

.987 

,998 

.998 

.  798 

-A.o 

IP.O 

.987 

.998 

.998 

.  798 

,987 

,996 

.998 

.798 

-1?.0 

IP.O 

.987 

1  .000 

1  .000 

.  798 

.987 

.998 

.998 

.798 

-18.0 

IP.O 

.987 

1  .000 

1  .000 

.798 

.987 

.998 

.998 

.798 

•PO.O 

IP.O 

.987 

1 ,000 

1 .000 

.  798 

0.000 

0,000 

0.000 

0.000 

-?a.o 

IP.O 

.987 

1  .000 

1  .000 

.798 

0.000 

0.000 

0,000 

0  .  ooo 

-?R.O 

IP^O 

.987 

1 .000 

1  .000 

.798 

0.000 

0.000 

0.000 

0.000 

IP^O 

.987 

.998 

.998 

.798 

0.000 

0.000 

0.000 

0.000 

-^8.0 

IP.O 

.987 

1 .000 

1  .000 

.798 

0.000 

O.OOO 

0,000 

0.000 

A-3 


-ao^o 

i;>^o 

0.000 

1  .000 

1 .000 

0.000 

0.000 

0.000 

0.000 

0.000 

•  au .  0 

1?.0 

0.000 

1  .000 

1  .000 

0.000 

.987 

UOOO 

1  .000 

.798 

•UR.O 

1?.0 

0.000 

1.000 

1  .000 

0.000 

0.000 

1  .000 

1.000 

0.000 

-s?.n 

1?^0 

0.000 

1  .000 

1  .000 

0,000 

0.000 

0.000 

0.000 

0,000 

1?,0 

0.000 

1.000 

1.000 

0.000 

0.000 

0^000 

0.000 

0.000 

•hO,  0 

1P.0 

0.000 

1  .000 

1  .000 

o.ono 

0,000 

1 .000 

1  .000 

0.000 

-ha.o 

IP^O 

O.ooo 

1  .000 

1  .000 

o.ooo 

0,000 

1  .000 

1  .000 

0,000 

0 

1P.0 

0.000 

1 .000 

1^000 

0.000 

0.000 

1  .000 

I  .000 

0.000 

-7P^n 

1P.0 

0.000 

1  .000 

1  .000 

0.000 

0.000 

1.000 

1  .000 

0.000 

-7^>.0 

IP.O 

0.000 

1  .000 

1,000 

0 .000 

o.ooo 

1  .000 

l.OOO 

0.000 

-RO.O 

IP.O 

o.ooo 

1  .000 

1  .000 

o.ooo 

o.ooo 

1.000 

1 .000 

0.000 

-flU.O 

1?.0 

0.000 

1.000 

1  .000 

0.000 

0.000 

1 .000 

1  .000 

0.000 

-HR.O 

1P.0 

0.000 

1 ,000 

1.000 

0.000 

0.00  0 

1^000 

1 .000 

0.000 

-9?.0 

iP.n 

0.000 

1  .000 

1  .000 

0.000 

0.000 

1^000 

r.aoo 

0,000 

•96 . 0 

IP^O 

0.000 

1  .000 

1.000 

0.000 

0.000 

l.OOO 

l.OOO 

0.000 

100. 0 

1P.0 

0.000 

1 .000 

1  .000 

0.000 

0.000 

1.000 

1.000 

0.000 

tO/J.O 

IP.O 

o.ooo 

1  .000 

1  .000 

0.000 

0.000 

1 .000 

l.OOO 

0.000 

1  0  R  .  0 

1P.0 

0.000 

1  .000 

1  .000 

0.000 

0.000 

1 .000 

1 .000 

0,000 

1  1?.0 

1P.0 

0.000 

1  .000 

1  .000 

0,000 

o.ono 

1^000 

1  .000 

0,000 

116,0 

1P.0 

0.000 

1 .000 

1  .000 

0,000 

0.000 

1  .000 
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